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Concentrated  sulfuric  acid  and  oleum  are  among  the  most  potentially 
hazardous  chemicals  routinely  transported  in  bulk  quantities  on  U.S.  and 
international  waterways.  Conceivably,  during  a  marine  mishap,  tons  of 
sulfuric  acid  could  abruptly  be  released  into  the  water,  and  the 
consequences  of  such  a  spill  accident  could  be  detrimental  to  man  and  the 
environment.  In  the  present  report,  several  acid  spill  accidents  are 
sunmarized.  It  is  shown  that  the  convective  mixing  of  of  concentrated 
sulfuric  acid  with  water  can  adequately  be  described  by  the  variation  of  the 
buoyancy  force  arising  from  changes  in  acid  concentration  and  released  heat 
of  dilution.  Mathematic  models  are  presented  for  the  cases  of  both 
instantaneous  and  continuous  releases.  For  oleum  spills  in  which  acid 
aerosol  formation  is  a  potential  safety  hazard,  a  conservative  estimate  of 
less  than  one-tenth  of  a  percent  is  obtained  for  the  amount  of  acid  which 
might  become  airborne  under  most  accident  conditions.  The  fraction  of 
airborne  acid,  however,  decreases  very  rapidly  with  increasing  release  depth 
below  water  surfaces.  The  acid  aerosols  exhibit  a  well-defined  log-normal 
particle-size  distribution  with  peak  diameter  varying  from  0.1  to  0.6  m(at 
70%  R.H.)  depending  upon  release  depths.  This  is  well  within  the  respirable 
particle  size  range. 


17.  Key  Vprds 

Sulfuric  Acid,  Marine  Spill, 
Acid  Aerosol 


II.  DiiiriNtiw  Statement 


Document  available  through  National 
Technical  Information  Service, 
Springfield,  Virginia  22161 


If.  ItMAI,  Clnti).  W  'W’l 

X.  Sawrlty  Claatit.  (•*  «fcia  *•*•) 

Unclassified 

Unclassified 

SULFURIC  ACID  SPILL  CHARACTERISTICS  UNDER 
MARITIME  ACCIDENT  CONDITIONS 


A  Final  Report 


prepared  for 


United  States  Coast  Guard 
Department  of  Transportation 


by 


I.N.  Tang,  U.T.  Unna. 

H.R.  Munkelwitz  and  G.S.  Smith 
Department  of  Energy  and  Environment 
Brookhaven  National  Laboratory 
Upton,  New  York  11973 


This  work  was  performed  under  the  auspices  of  the  United  States  Coast  Guard 
under  Contract  No.  MIPR-2-70099-7-75117A. 


v 


CONTENTS 


Introduction  . 1 

Sulfuric  Acid  Spill  Scenarios  .  2 

Class  I  Accidents:  Convective  Mixing  of  Acid  with  Water  .  4 

1.  Instantaneous  Releases  .  4 

(1)  Theoretical  Considerations  . . . . .  4 

(ii)  Experimental  Results  and  Discussions  .  8 

2.  Continuous  Releases  . .  11 

(1)  Theoretical  Considerations  .  11 

(ii)  Experimental  Results  and  Discussions  .  18 

Class  II  Accidents:  Acid  Aerosol  Formation  .  21 

Conclusion  . 26 

References  .  27 

Appendix  1.  Computer  Program  for  Instantaneous  Releases  .  29 

Appendix  2.  Photographs  of  the  BNL  Acid  Spill  Facility  .  30 

Appendix  3.  Evaluation  of  the  Integrals  in  Eq.  (15)  to  (18)  .  31 

Appendix  4.  Computer  Program  for  Continuous  Releases  .  32 

Appendix  5.  A  Non-linear  Iterative  Algorithm  for  Inverting 

Diffusion  Battery  Data  . 33 


ii 

* 

-  —  -  *-  - -  — 


INTRODUCTION 


Concentrated  sul£urlc  acid  and  oleum  are  among  the  most  potentially  haz¬ 
ardous  chemicals  routinely  transported  in  bulk  quantities  on  U.S.  and  interna¬ 
tional  waterways.  Although  these  chemicals  are  normally  contained  in  cargo 
compartments  sealed  against  gas  and  water  leakage,  the  rupture  of  one  or  more 
of  these  compartments  during  a  marine  accident  could  conceivably  bring  about 
the  abrupt  release  of  tons  of  acid  into  water.  The  large  amount  of  heat  gen¬ 
erated  from  acid-water  mixing  could  cause  both  water  and  the  absorbed  sulfur 
trioxide  in  oleum  to  vaporize,  thus  forming  an  acid  mist  in  the  atmosphere. 

This  acid  mist  would  pose  an  immediate  danger  to  anyone  directly  involved  in 
the  accident  and,  under  adverse  meteorological  conditions,  even  threaten  the 
safety  of  the  nearby  public  as  well.  Furthermore,  depending  upon  the  magnitude 
and  location  of  a  specific  spill,  the  water  temperature  and  acidity  could  be 
significantly  altered  long  enough  to  destroy  marine  life  in  the  release 
vicinity. 

There  have  been,  in  the  past,  several  documented  marine  accidents  involv¬ 
ing  sulfuric  acid  cargo  vessels.  The  most  recent  acid  spill  accident,  known  as 
the  case  of  the  KT  Big  Mama  [1],  occurred  at  dawn  on  August  18,  1976.  A  tank 
barge  loaded  with  approximately  1,060  tons  of  20%  oleum  capsized  while  being 
towed  in  the  lower  Chesapeake  Bay.  All  of  the  cargo  was  lost  into  the  Bay  dur¬ 
ing  a  time  period  of  30  to  60  min.  The  towboat  personnel  later  recalled  that 
they  observed  a  tremendous  reaction  under  the  capsized  barge,  which  produced 
great  quantities  of  steam  and  vapor.  In  fact,  it  was  noted  that  the  barge 
itself  was  lifted  4  feet  upward  in  the  water  for  an  estimated  30  minutes  before 
sinking.  Fortunately,  the  accident  took  place  in  a  remote  deep-water  region  of 
the  Bay  and  no  apparent  damage  was  reported  except  for  a  rather  extensive  but 


localized  fish  kill  sighted  by  the  first  Coast  Guard  boatcrew  arriving  on  the 
scene. 

In  view  of  the  potential  hazard  associated  with  a  probable  sulfuric  acid 
spill,  the  United  States  Coast  Guard  (USCG)  initiated  a  safety  analysis  program 
and  contracted  the  Brookhaven  National  Laboratory  to  carry  out  laboratory 
experiments  designed  specifically  for  investigating  (a)  the  mechanisms  by  which 
sulfuric  acid  aerosols  may  be  formed  from  oleum  spills  on  water,  and  (b)  the 
mixing  processes  which  take  place  in  water  within  the  proximity  of  a  release 
source.  The  results  of  the  study  are  reported  here,  together  with  mathematical 
models  developed  during  the  course  of  investigation  for  predicting  acid  disper¬ 
sion  in  water.  The  quantitative  information  presented  here  provides  a  data 
base  for  USCG  to  set  up  emergency  response  procedures  and/or  make  environmental 
impact  statements  in  case  of  a  sulfuric  acid  spill  accident. 

SULFURIC  ACID  SPILL  SCENARIOS 

For  simplicity  in  describing  the  acid  spill  scenarios,  the  following  two 
broad  classes  of  accidents  are  considered:  Class  I,  accidents  during  which  no 
appreciable  amount  of  acid  aerosols  are  produced  and  Class  II,  accidents  in 
which  the  major  safety  hazards  are  associated  with  acid  aerosol  formation.  In 
all  cases,  of  course,  a  large  amount  of  heat  is  liberated  from  acid-water 
mixing,  especially  during  the  very  early  stage  of  the  spill. 

Thus,  the  Class  I  accidents  Involve  only  the  convective  mixing  of  acid 
with  water,  resulting  from  either  a  quick  discharge  of  concentrated  sulfuric 
acid  into  deep  water,  such  as  might  be  the  case  of  a  capsized  barge  with  all 
hatches  blown  open,  or  a  gradual  leak  through  a  hole  of  finite  size  developed, 
say,  in  a  minor  collision  mishap.  The  third  possibility  is  the  spill  of  an 
oleum  cargo  well  underwater  as  in  the  case  of  the  MT  Big  Mama  described  above. 
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In  all  cases,  Che  released  acid  mass  tends  to  sink  quickly  in  water  because  of 
its  high  density.  Vigorous  turbulent  mixing  takes  place  over  the  advancing 
front,  resulting  in  entrainment  of  the  ambient  fluid.  The  acid-containing 
region  Initially  grows  extremely  fast  by  this  entrainment  process.  As  the 
volume  becomes  larger  and  the  density  difference  smaller,  the  advancing  speed 

of  the  acid  mass  as  a  whole  slows  down  accordingly. 

c> 

The  Class  II  accidents  include  spills  of  oleum  over  water  surfaces  or  into 
shallow  water.  The  heat  generated  in  the  reaction  dissipates  via  two  different 
processes:  raising  the  temperature  of  the  mixing  mass  and  providing  latent 
heat  for  sulfur  trioxide  and  water  vaporization.  The  evolution  of  SO3  is 
expected  to  occur  only  very  briefly  until  enough  water  is  entrained  to  dilute 
the  oleum  down  to  essentially  100%  H2SO4.  At  this  point,  the  "free"  SO3  is  all 
consumed  and  vaporization  of  the  H2SO4  is  unlikely  because  of  its  extremely  low 
volatility.  The  initially  vaporized  SO3  reacts  quickly  with  water  vapor  to 
form  sulfuric  acid,  which  will  then  nucleate  and  grow  into  an  acid  mist.  The 
acid  mist  will  then  be  carried  aloft  either  by  the  hot  buoyant  air  or  dispersed 
by  the  prevailing  wind. 

It  should  be  recognized,  however,  that  accident  conditions  are  neither 
predictable  nor  controllable,  and  that  the  spill  scenarios  described  above  are 
nothing  more  than  hypothetical  cases  in  which  each  stage  of  an  accident  may  be 
clearly  defined.  The  purpose  is  to  delineate  the  complex  phenomena  with  sev¬ 
eral  Individual  processes  which  are  reasonably  decoupled  for  the  convenience  of 
laboratory  investigation  over  a  variety  of  conditions.  Since  no  accident  may 
be  regarded  as  being  typical,  the  emphases  are  then  placed  on  those  extreme 
conditions  whose  outcomes  provide  a  conservative  estimate  for  most  other  cases. 
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Class  I  Accidents;  Convective  Mixing  of  Acid  with  Water 
1.  Instantaneous  Releases 

(1)  Theoretical  Considerations 

The  notion  of  an  acid  mass,  suddenly  released  in  water  as  the  result  of  a 
spill,  has  some  features  in  common  with  that  of  a  thermal  rising  through  stag¬ 
nant  air.  In  both  cases,  the  essential  driving  force  originates  from  the 
presence  of  a  buoyancy  source  due  to  the  density  difference  between  the  moving 
fluid  and  the  ambient  fluid.  The  physics  of  thermals  has  been  studied  rather 
extensively  since  the  1950s  by  many  Investigators  including  among  others, 
Batchelor  [2],  Scorer  [3],  Turner  [4],  Csanady  [5],  and  Tsang  [6].  In  the 
present  study,  the  simple  model  of  Morton  et  al.  [7]  is  adapted  for  the  case  of 
acid-water  mixing.  An  important  modification  of  the  model  is  the  incorporation 
of  an  energy  equation  to  account  for  the  simultaneous  dilution  and  heat  genera¬ 
tion.  The  released  heat  tends  to  induce  further  turbulence  by  creating  temper¬ 
ature  inhomogeneities  in  addition  to  concentration  inhomogeneities  arising  from 
entrainment  processes. 

Consider  a  lump  of  concentrated  sulfuric  acid  discharged  instantaneously 
into  sti 11  water.  A  sequence  of  the  event  is  shown  photographically  in  Fig. 

1.  As  seen  in  these  pictures,  the  boundary  of  the  acid  mass  is  highly  irregu¬ 
lar  due  to  the  turbulent  motion  both  within  and  adjacent  to  the  settling  mass. 
However,  an  equivalent  sphere,  with  radius  b’(t')  which  is  a  function  of  time, 
t',  may  be  defined  at  any  instant  to  envelop  the  bulk  of  the  acid  mass.  Let 
x'(t')  represent  the  position  of  the  center  of  the  sphere  relative  to  the  point 
of  release  and  u'(t')  the  mean  velocity  of  the  center.  Hereafter,  the  primed 
notation  refers  to  quantities  expressed  in  the  c.g.s.  units,  whereas  the 
unprimed  notation  refers  to  dimensionless  quantities  to  be  defined  as  they 
appear. 
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The  key  assumption  of  the  model  Is  that  the  rate  of  entrainment  is  propor¬ 


tional  to  u'(t').  The  governing  equations  of  motion  may  then  be  written  as 
follows : 

Mass  balance: 

_fL_  irb,3p')  ■  Aitb'^au'p*  (1) 

dt'  3  ' 

Momentum  balance: 

_L  (l_*b'Vu')  -  1  irb,3g(p*  -p  ')  (2) 

dt'  3  3  W 


Energy  balance: 


Acid  balance: 


(i  ^b'3p'y) 


0 


Center-of-mass  velocity: 

dx'  , 

- -  ■  u 

dt’ 


(3) 


(4) 


(5) 


where  p'  and  p^  are  the  density  of  the  acid  sphere  and  ambient  fluid, 
respectively,  y(t’)  the  average  H2SO4  mass  fraction  of  the  solution  in  the 
sphere,  (V  the  specific  heat  at  concentration  y,  AH'  the  heat  of  dilution  per 
unit  mass  of  the  solution  at  concentration  y,  T'(t')  the  excess  temperature  of 
the  sphere  over  the  ambient  fluid,  and  g  the  gravitational  constant.  Eq .  (1) 
expresses  the  conservation  of  total  mass  (or  equivalently,  the  mass  balance)  as 
the  acid  sphere  expands.  Eq .  (2)  is  the  conservation  of  momentum.  Eq .  (3)  is 
the  conservation  of  energy.  Eq.  (4)  is  the  conservation  of  the  total  H2SO4 
content.  Eq.  (5)  is  simply  the  definition  of  the  velocity,  u',  at  the  center 
of  the  acid  sphere. 
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It  is  a  common  practice  in  hydrodynamic  formulations  to  make  free  use  of 
the  so-called  Boussinesq  approximation,  namely,  the  density  difference  between 
the  bouyant  fluid  and  the  ambient  fluid  is  assumed  negligible  except  in  the 
buoyancy  term.  Since  the  initial  density  of  the  96%  sulfuric  acid  is  1.836 
gm/cm^  at  room  temperature,  it  can  hardly  be  regarded  as  being  close  to  the 
ambient  water  density.  Thus,  in  the  present  case  the  Boussinesq  approximation 
is  not  justifiable  especially  in  the  early  stage  of  the  motion.  We  have  there¬ 
fore  retained  the  density  of  the  acid  mass  on  the  left  hand  sides  of  Eqs .  (1)  - 
(A).  The  resulting  governing  equations  are  more  complicated  than  would  be  the 
case  if  the  Boussinesq  approximation  had  been  applied,  because  '  is  a  function 
of  y  and  T'  and  hence  Eqs.  (1),  (3),  and  (4)  are  coupled.  Clearly,  the  effect 
of  finite  density  difference  between  the  acid  mass  and  ambient  fluid  is  partic¬ 
ularly  important  in  the  momentum  and  energy  equations. 

For  the  case  where  the  acid  mass  is  released  with  zero  initial  momentum, 
the  initial  conditions  are, 

x'  -  o,  u'  -  o,  T*  -  T'0>  y  -  y0,  b’  =  b^  (6) 

where  b^  and  y0  are  the  equivalent  radius  and  H2SO4  mass  fraction  of  the  start¬ 
ing  sample,  respectively.  Since  in  all  experiments  the  concentrated  sulphuric 
acid  is  usually  warmer  than  the  ambient  water  temperature,  T'0  is  taken  to  be 
5°C  throughout  the  study. 

The  specific  heat  data  as  a  function  of  acid  composition  at  20°C  were 
taken  from  Perry  and  Chilton  [8].  Since  literature  data  [9-11]  indicate  only  a 
very  small  temperature  dependence  for  the  specific  heat  of  sulfuric  acid,  the 
rather  comprehensive  room  temperature  data  of  Reference  [8]  were  therefore  used 


in  the  calculation: 


cp(y)  = 


1  -  0.791y 
0.603  -  0.192y 
1.109  -  0. 791y 


0 <  y  _<  0.6625 
0.6625  <  y  <  0.8448 
0.8448  <  y 


The  density  and  heat  of  dilution  data  were  computed  from  the  empirical 
relationships  given  by  Azarniouch  et  al. (12)  as  follows: 


m)f(T’)-i 


-20°C  <  T*  <  80°C, 

where  <p  (£)  =  1.0  +  3.0943  £  -  4.0834£2  +  1.9032?3 


5.444  -  4.444y 


f(T')  -  1.0  -  0.000536  T’ 

and 

(i^l)  =  -  (A  -  By  +  Cy2),  cal/g 

where  A  «  492.7,  B  «  987.9  and  C  =  558.3.  Again,  the  effect  of  temperature  on 
heat  of  dilution  is  neglected  here. 

To  facilitate  numerical  integration  of  Eq.  (1)  through  (5),  the  variables 
x',  b',  u',  y  and  Tf  were  first  grouped  as  follows: 

V'  -  b'3,  M'  -  b'3u',  W'  -  Cp'b’V,  N'  -  b’3y. 


These  were  then  converted  into  the  dimensionless  groups: 


V  -  v0v, 


M 


M, 


1/2  »7/6, 
vo 

1/3 


t’ 


W  -  (AVoy0)W, 

-  (e"1/2  v;1/6)t 
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N’  -  (y0Vo)N,  x-  -  v; 


X 


where  \M  -  b^.  By  changing  variables  accordingly,  incorporating  the  empirical 
expressions  for  C^,  p*  and  3AH'/3y,  and  rearranging,  Eq.  (1)  through  (5)  were 
transformed  into  dimensionless  forms  suitable  for  simultaneous  integration. 

The  corresponding  dimensionless  differential  equations  are  as  follows: 


dV 

dt 


dM 

dt 

dW 

dt 

dN 

dt 


3c,  M  ,  .  -  I.  ($f) 

1/3  $f  dt 


4>fV 


vl£-i 

<|>f  4>f  dt 

d  N  N  N  2  u  d 

-  (I)  [l-Xjil  +  x2<->  1  ■  ^  4-  (♦*> 

dt  V  1  V  V  <j>  f  dt 


(<ff) 

4»f  dt 


dx  m  M 
dt  V 


(7) 

(8) 

(9) 

(10) 

(ID 


2 

where  Xj  «  By0/A  and  X2  ■  Cy0/A. 

The  numerical  integration  of  Eqs.  (7)  through  (11)  were  carried  out  on  a 
CDC  6600,  using  the  fourth-order  Runge-Kutta  numerical  algorithm.  The  initial 
conditions  are: 


V  -  1,  M  -  0,  W  -  CP(yo)T°  ,  N  -  1,  x  -  0 

Ay0 

A  print-out  of  the  computer  program  for  instantaneous  release  calculations 
is  given  in  Appendix  1 . 

(ii)  Experimental  Results 

Experiments  designed  to  study  both  acid-water  mixing  and  acid  aerosol  for¬ 
mation  under  simulated  accident  conditions  were  carried  out  in  the  Acid  Spill 
Facility,  which  is  shown  schematically  in  Fig.  2.  Basically,  the  apparatus 
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consists  of  a  4'  x  4*  x  8*  stainless  steel  tank,  equipped  with  various  acid 
release  devices,  temperature  probes,  and  an  aerosol  sampling  system.  One  side 
of  the  tank  was  constructed  entirely  of  a  transparent  Lexon  panel  for  photo¬ 
graphic  purposes.  A  smaller  2'  x  2*  x  4'  lucite  tank  was  used  initially  in  the 
study.  Some  photographs  of  the  Acid  Spill  Facility  are  given  in  Appendix  2. 

Several  different  types  of  acid  release  devices  were  used  in  the  experi¬ 
ments.  In  the  first  series  of  experiments,  hemispherical  glass  cups  (shown  in 
Fig.  3),  designed  to  pivot  about  a  horizontal  axis  fixed  along  its  diameter, 
were  used  to  study  the  instantaneous  release  of  cone,  sulfuric  acid  in  water. 
Attached  to  the  cup  was  a  glass  lever  with  a  float  on  one  end,  which  when 
placed  in  water  together  with  the  cup,  served  to  stabilize  the  cup  before 
inverting  it.  In  a  typical  experiment,  the  cup  was  filled  with  a  known  amount 
of  acid  to  the  level  of  the  ambient  water.  It  was  then  rapidly  inverted  by  a 
specially  designed  mechanical  device  to  release  the  acid  mass  without  any  ver¬ 
tical  Impulse.  After  dumping,  the  cup  quickly  reverted  to  its  original  posi¬ 
tion,  thereby  eliminating  the  problem  of  post -dripping.  By  properly  illuminat¬ 
ing  the  water  tank  with  an  Intense  tungsten  lamp,  the  motion  of  the  acid  mass 
was  readily  observable  because  of  the  refractive  index  difference  between  the 
two  fluids.  Shadowgraphical  pictures  of  the  settling  acid  mass  were  taken 
against  a  Scotchcal  film  background  at  various  times  after  release  by  a  camera 
equipped  with  an  automatic  shutter  release  and  film  advance  device.  Fig.  4 
shows  a  series  of  shadowgraphs  taken  at  various  stages  of  mixing  for  a  500  g 
acid  dump.  Several  selected  blow-ups,  shown  earlier  in  Fig.  1,  clearly  indi¬ 
cate  that  the  acid  mass  retains  a  reasonably  spherical  shape  throughout  the 
course  of  its  motion,  thus  somewhat  lending  support  to  the  equivalent  sphere 
assumption  used  in  the  mathematical  model. 
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In  experiments  carried  out  to  study  the  dispersion  of  concentrated 
sulfuric  acid  (96X  by  wt.)  in  water,  four  different  initial  sample  sizes 
spanning  two  orders  of  magnitude  in  weight  were  employed.  These  were  4.6,  22, 
130  and  500  gm.  Several  runs  were  made  for  each  sample  size  and  the  results 
from  all  experiments  are  shown  in  Fig.  5,  in  which  the  dimensionless  frontal 
distance  of  the  settling  acid  mass  given  by  (x*  +  b'J/bg  is  plotted  vs.  the 
dimensionless  time,  t ' V  g/bg,  for  any  given  initial  sample  size  bg.  It  is 
evident  from  Fig.  5  that  all  data  points  gather  around  a  single  curve  regard¬ 
less  of  the  initial  sample  size  used.  The  solid  curve  in  Fig.  5  represents  the 
"best  fit"  results  calculated  from  the  mathematical  model  with  a  »  0.21. 

Morton  et  al.  (7)  have  conducted  similar  experiments  by  suddenly  releasing  a 
puff  of  light  fluid  in  a  stably  stratified  salt  solution  which  is  denser  than 
the  rising  fluid.  They  found  that  the  values  of  a  lie  between  the  limits  0.27 
and  0.34,  which  are  certainly  larger  than  the  0.21  value  obtained  in  the 
present  study  for  the  turbulent  mixing  of  sulfuric  acid  and  water.  The  differ¬ 
ence  can  perhaps  be  attributed  to  the  heat  of  dilution  that  is  accompanying  the 
mixing  process  in  the  present  case  and  to  the  non-Boussinesq  behaviour. 

It  is  of  interest  to  know  how  the  average  temperature  of  the  acid  mass 
varies  with  time.  Attempts  were  made  to  measure  local  temperatures  with  fast- 
response  (0.03  sec  time  resolution)  miniature  thermocouples  placed  at  different 
locations  through  which  the  settling  acid  mass  was  expected  to  pass.  Due  to 
the  turbulent  nature  of  the  dispersing  acid  mass,  however,  large  local  tempera¬ 
ture  fluctuations  were  observed  from  run  to  run,  and  it  was  extremely  difficult 
to  Interpret  the  measured  temperature  data.  Instead,  the  characteristic  excess 
temperatures,  T' ,  as  a  function  of  dimensionless  time  were  calculated  using  the 
mathematical  model,  and  Fig.  6  shows  the  results  obtained  with  a  -  0.21.  It 
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is  seen  that  the  excess  temperature  rises  to  a  maximum  very  shortly  after 
release  and  then  decays  rather  rapidly  with  time.  The  calculated  temperature 
variation  with  time  was  qualitatively  confirmed  with  measurements.  Also  shown 
in  Fig.  6  is  the  change  in  acid  mass  fraction  as  a  function  of  dimensionless 
time.  For  example,  a  100-fold  dilution  is  expected  to  occur  in  3.8  sec  for  an 
initial  acid  size  of  500  gm,  and  8.3  sec  for  a  50  kg  release. 

2.  Continuous  Releases 

(1)  Theoretical  Considerations 

The  dispersion  of  sulphuric  acid  from  a  gradual  leak  through  a  finite  size 
hole  resulting  from  a  minor  vessel  collision  may  be  simply  described  in  terms 
of  existing  models  for  an  axisymmetric  maintained  plume  in  a  uniform  quiescent 
fluid.  A  large  volume  of  work  exists  on  modeling  buoyant  plumes  and  our 
present  state  of  knowledge  appears  to  be  adequate  in  describing  their  behavior 
(for  good  reviews  see  for  examples  {13-15]).  However,  because  of  the  high 
pressure  head  in  the  acid  container  of  the  vessel  (as  much  as  1  kiloton),  the 
acid  will  emerge  from  the  hole  with  large  momentum  flux.  Such  a  scenario  is 
characterized  by  some  researchers  as  a  forced  plume  [16]  or  as  a  buoyant  jet  by 
others  {17].  The  present  model  will  somewhat  resemble  Morton's  formulation 
[16].  The  general  characteristics  of  the  model  are  depicted  schematically  in 
the  following  sketch: 
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Underlying  this  model  are  the  following  assumptions: 

(1)  steady-state  condition  prevails  so  that  all  the  physical  variables 
are  functions  of  x'  only. 

(2)  ambient  fluid  is  uniform  and  quiescent. 

(3)  the  plume  is  axisymmetric  and  turbulence  is  fully  developed  except 
very  close  to  the  source. 

(4)  the  flow  is  self-similar.  Accordingly,  the  velocity  profile  is 
assumed  to  be  Gaussian, 

^  _ri2/ui2 

u’(x\  r')  -  u(x’)e  r  /D  (12) 

where  b'  -  b'(x')  may  be  regarded  as  the  characteristic  radius  of  the 
plume  at  x'  and  u(x')  is  the  axial  velocity. 

(5)  Experimental  data  indicate  that  the  temperature  profile  is  also 
Gaussian  but  with  a  larger  lateral  spread  than  that  of  velocity, 

~  -r’2/x2b'2 

T'(x' ,r’)  -  T(x’)e  /A  D  (13) 

where  T'  stands  for  excess  temperature  over  that  of  the  ambient  fluid 
and  O  1  Is  an  empirical  constant.  Similarly,  the  acid  concentra¬ 
tion  is  expected  to  be  higher  in  the  core  of  the  plume  than  at  the 
"outer"  edge.  Moreover,  the  temperature  rise  is  caused  by  the  heat 
generated  due  to  acid  dilution.  Therefore,  the  acid  mass  fraction 
distribution  is  assumed  to  be  Gaussian  with  the  same  lateral  spread 
as  that  of  the  temperature  profile. 
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~  -r'2/x2b'2 

y’Cx'.r')  -  y(x’)e  ' 


(6)  For  computational  conveniences,  here  the  Boussinesq  approximation  is 
retained,  i.e.  the  density  of  the  plume  is  assumed  constant  and  equal 
to  that  of  the  ambient  fluid  except  in  the  buoyancy  term 


-r ' 2/\2 wi 2 

gAp’(x',r‘)  -  gp;P(x')e  M  b 


(Ha) 


where  P(x’)  -  ^(y)f(T)  -  1 


(7)  Although  the  heat  capacity  of  the  acid  is  a  function  of  y*  and 

therefore  would  have  a  lateral  spread  too,  it  is  assumed  to  be  given 
by  the  mass  fraction  y  at  the  axis  of  the  plume 


Cp(x')  -  Cp(y(x' )) 

(8)  The  rate  of  entrainment  of  ambient  fluid  at  the  "edge”  of  the  plume 
at  any  level  is  proportional  to  the  axial  velocity  at  that  level. 

The  proportionality  constant,  a,  is  Independent  of  x' . 

With  these  assumptions  the  following  four  conservation  equations  can  be  written 


volume:  - —  (  JYirr'u’dr*  J  -  2itbr 

dx'  \o  / 


momentum: 


y  CD  V  CD 

i — (  fiirr'u'Zfi'dr')  •  r2«r,gAp'dr' 
dx'  \o  /  o 
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(b'2uC  T) 
dx'  y 


-  -<x2+i)b,2u  [?(?)£  +  S(?)il 

L  dx  b'  dx'J 


_  (b,2uy) 


The  details  of  the  derivation  are  given  In  Appendix  3. 

To  transform  Eqs.(19)-(22)  into  non-dimensional  forms,  we  need  to  choose 
the  correct  scales  for  the  physical  variables  that  govern  the  behavior  of  the 
plume.  At  the  source  the  momentum  force  is  expected  to  dominate  over  buoyancy 
force  and  therefore  the  momentum  flux  M0  ■  b,2u2  at  the  source  is  of  primary 
Importance  in  setting  the  scales.  Since  b'  is  a  statistical  variable,  M0  must 
be  properly  related  to  measurable  quantities, 

M0  ■  (b,2u2)x,_Q  ”  (volumetric  flow  rate)  x  (mean  velocity) 


mass  flow  rate 
density  of  fluid 


x  (mean  velocity) 
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For  a  source  in  the  form  of  a  circular  orifice  with  radius  R, 


mean  velocity  » 


mass  flow  rate 
(density)irR2 


Therefore, 


mass  flow 
density  x 


rate  |  ■ 

~  J 


Then  the  following  transformation  is  chosen 


b,2U  -  (M05/6g_1/3)V 
b,2u2  -  M0M 

b'^CpT  -  (YM05/6g“1/3Cp*)W 
b,2uy  -  (y0Mo5/,68  1^3)N 
x’  -  (M^V'3* 
b’  -  (Mo1/3g"1/3)b 

Where  y  -  Ay0/Cp*.  Then,  in  non-dimensional  form,  Eqs.  (19)-(22)  become 


£!  -  2oM1/2 


—  -  2X2l-<4)f-l) 
dx  M 


(A2*!)*  r(Q-2S)^  +  ^£ll 
dx  V  [  dx  M  dx  J 


0 


where 


It  must  be  emphasized  that  Eqs.  ( 24 ) — ( 27 )  are  good  approximations  only  for 
a  fully  developed  free  turbulent  plume  at  a  distance  far  from  the  source. 

Close  to  the  source  the  situation  is  more  complicated  and  the  nature  of  the 
flow  1*  evolving  very  rapidly,  due  to  Instability,  from  a  laminar  Jet  immedi¬ 
ately;  after  emergence  from  the  source  to  a  fully-developed  self-preserving  tur- 
buler'  jet  at  some  distance  x^'.  To  construct  a  model  which  is  valid  for  the 
entire  plume  is  difficult  and  besides,  given  the  complexity  of  an  actual  acid 
spill,  it  is  probably  not  worth  the  effort.  Instead,  the  following  approach  is 
adopted.  The  flow  is  divided  into  two  regimes:  0_<  x’  _<  X]/,  the  flow  is 
essentially  a  pure  jet,  where  momentum  forces  dominate  and  for  x'  >  Xj',  it  is 
essentially  a  buoyant  plume,  where  the  buoyancy  force  becomes  dominant.  The 
governing  equations  ( 24)-( 27 )  are  assumed  to  be  valid  for  both  regimes  with 
different  values  of  the  entrainment  constant  a.  The  point  of  transition 
between  these  two  regimes,  xj',  is  somewhat  uncertain.  Tennekes  and  Lumley 
[18]  indicated  that  the  mean  velocity  profile  appears  to  be  self -preserving 
beyond  about  eight  orifice  diameters,  although  turbulence  quantities  are  still 
evolving  at  40  orifice  diameters.  Therefore,  given  the  simplification  of  the 
present  model,  X|'  will  be  adjusted  between  10  and  15  orifice  diameters  to 
yield  the  best  fit  with  measurements. 


Experimental  data  presented  by  Ricou  and  Spalding  [19]  Indicated  that  the 
entrainment  rate  for  a  buoyant  jet  is  somewhat  larger  than  that  for  a  pure  jet 
with  the  same  momentum*  Comparison  of  boundary- layer  theory  with  experiments 
showed  that  for  an  axisymmetric  pure  jet  the  volumetric  flow  rate  at  a  distance 
x'  from  the  orifice  may  be  given  by  [20] 

volumetric  flow  rate  =  0.404-^|Mo  x’  (28) 

Thus ,  the  gradient  of  the  volumetric  flow  rate  is  constant .  This  fact  may  be 
used  to  define  an  equivalent  entrainment  constant  for  a  pure  jet  based  on 
Taylor’s  entrainment  assumption  [7], 

_  (volumetric  flow  rate)  *  2Trab'u  (29) 

dx’ 


where  (b'u)2  is  a  measure  of  the  average  momentum  flux.  If  the  pure  jet  regime 
is  small  so  that  the  buoyancy  force  has  a  negligible  effect,  b’u  is  approxi¬ 
mately  constant  by  virtue  of  the  conservation  of  momentum.  Then,  comparison 
between  Eqs.  (28)  and  (29)  implies  that 


a 


0.404 

~U 


0.065 


(30) 


Direct  measurements  of  entrainment  into  a  pure  jet  yielded,  however,  a  smaller 
value  of  0.057  for  o  [17,21].  Thus,  the  value  of  a  will  be  adjusted  in  our 
model  between  0.057  and  0.065  to  give  a  best  fit  with  measurements.  For  x’  > 

’ ,  i.e.  in  the  buoyancy  regime,  existing  data  showed  that  a  lies  between  0.08 
and  0.12  [13,16,17,21];  the  lower  value  corresponds  to  Gaussian  profiles  while 
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the  higher  value  to  "top-hat"  profiles.  In  the  present  model,  we  again  adjust 
a  between  this  range  to  yield  best  fit  with  experimental  data. 

It  can  be  readily  shown  from  Eqn.  (24)-(27)  that  as  x  — >°° 


which  agrees  with  known  results  deduced  from  pure  plume  models  [7,22].  This  is 

intuitively  expected,  because  far  from  the  source  the  acid  in  the  plume  has 

been  so  diluted  that  the  flow  approaches  the  conditions  of  a  pure  plume.  Thus, 

at  large  x,  the  plume  tends  to  approach  the  shape  of  a  circular  cone  subtending 

6ct 

an  angle  of  2  tan“^(— — )•  The  measurement  of  this  angle  offers  a  simple  way  to 
determine  the  entrainment  constant  in  the  buoyant  regime  (x1  >  xi'). 

(ii)  Experimental  Results 

A  forced  plume  of  96%  sulphuric  acid  was  produced  from  a  release  device 
constructed  simply  by  fusing  an  ordinary  laboratory  burett  to  the  bottom  of  a 
glass  trough.  The  tip  of  the  burett  was  placed  just  above  the  water  surface  in 
the  acid  spill  tank.  The  glass  trough  is  made  as  large  as  possible  so  that 
during  the  course  of  the  experiment  the  head  of  the  acid  in  the  trough  was 
fairly  constant  thereby  maintaining  a  steady  rate  of  acid  feed.  This  is  essen¬ 
tial  to  assure  steady-state  condition.  The  mass  flow  rate  was  measured  by 
weighing  the  amount  of  acid  released  in  a  certain  period  of  time  and  was  found 
to  be  about  22.2  gm/sec. 

High  speed  photographs  were  taken  at  various  Instants  after  the  acid 
release,  again  using  shadowgraphic  techniques.  Figs.  7a  and  7b  show  a  typical 
set  of  pictures  which  were  taken  at  a  rate  of  1  frame  per  second.  It  can  be 
seen  that  at  about  10  seconds  after  release,  steady-state  was  reached.  The 


acid  feed  was  turned  off  as  soon  as  the  acid  front  reached  the  bottom  of  the 
tank.  Fig.  8,  9,  and  10  are  enlargements  of  the  plume  of  the  same  run  at  13, 

21  and  29  seconds,  respectively,  after  the  release.  A  comparison  between  Fig. 

9  and  10  indicates  that  the  shape  and  the  size  of  the  plume  were  pretty  stable 
with  the  bulk  of  the  acid  flow  more  or  less  confined  within  a  circular  cone. 

Temperature  distribution  across  a  horizontal  section  of  the  plume  was 
measured  by  placing  an  array  of  four  thermocouples,  equal  distances  apart,  so 
that  the  first  thermocouple  was  vertically  below  the  release  nozzle,  i.e.  on 
the  axis  of  the  plume,  and  the  fourth  thermocouple  more  or  less  at  the  "edge" 
of  the  plume.  Fig.  11-1  shows  the  measurements  with  four  such  thermocouples 
placed  at  a  distance  two  inches  below  the  nozzle,  showing  a  distinct  tempera¬ 
ture  profile.  Fig.  11-2  indicates  that  the  temperature  profile  at  nine  inches 
from  the  nozzle  was  almost  flat.  Only  measured  plume  axis  temperature,  T' ,  as 
a  function  of  distance,  x',  from  the  nozzle  can  be  directly  compared  with  theo¬ 
retical  predictions.  The  experimental  data  are  summarized  in  Table  1. 


xf 

(inches) 

X 

T’ 

<°c> 

T  =  T’  -  Tw' 
(V  *  18°C) 

2 

3.96 

121.0 

103.0 

3 

5.94 

/2.0 

54.0 

4 

7.92 

62.0 

44.0 

5 

9.90 

56.0 

34.0 

6 

11.80 

46.0 

28.0 

7 

xi.86 

38.0 

20.0 

8 

15.84 

35.0 

17.0 

9 

17.82 

-  - 

33.0 

15.0 

Table  1.  Measured  plume  axis  temperatures  at  various 
distances  from  the  release  nozzle.  Column  2 
lists  the  dimensionless  distances  converted 
from  Column  1.  Column  4  contains  the  excess 
temperature  derived  from  Column  3. 
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With  a  measured  mass  flow  rate  of  22*2  g/sec,  Mq  was  calculated  to  be 
about  2069.5  cmVsec^,  which  defines  a  characteristic  length  scale  of  1.283  cm 
for  the  model.  From  Fig.  8,  9,  and  10,  the  angle  of  the  cone  defined  by  that 
part  of  the  plume  far  away  from  the  source  was  estimated  to  be  about  15°  ± 

0.5°,  which  would,  according  to  Eq.  31,  give  a  value  of  0.11  ±  0.004  for  the 
entrainment  constant  a  in  the  buoyant  plume  regime. 

The  proper  initial  conditions  for  the  model  are  at  x  *  0 

V  -  1 
M  -  1 
W  »  5/y 
N  -  1 

Here  T  *  5,  since  the  initial  acid  temperature  was  always  about  5°C  warmer  than 
the  water  temperature.  Since  Eq.  (27)  implies 

N  »  constant  •  1  ,  for  all  x  , 

the  governing  equations  reduce  to  three,  (24)-(26).  They  were  solved  simulta¬ 
neously  by  the  fourth  order  Runge-Kutta  algorithm  with  a  step  size  of  0.1  and 

a  -  ai  0  <_  x  £  x\ 

a  -  0.11  x  >  x^ 

A  print-out  of  the  computer  program  is  given  in  Appendix  4. 

The  values  of  aj  and  x^  were  adjusted  to  produce  a  best  fit  with  the 
temperature  measurements.  It  was  found  that  -  0.065  and  xj^  %  3.0  (which 
corresponds  to  about  13  orifice  diameters  from  the  source)  were  the  optimal 
figures.  In  Fig.  12,  the  theoretical  axial  excess  temperature  and  axial  acid 
mass  fraction  were  plotted  against  x  on  the  same  graph.  The  overall  agreement 
between  theoretical  prediction  and  experimental  data  for  axial  excess 
temperature  is  excellent. 
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With  Che  values  for  Che  entrainment  consCanc  In  boCh  regimes  escablished 
by  experimental  measurements,  Che  axial  acid  mass  fracClon  (as  shown  in  Fig. 

12)  was  predicted  Co  drop  below  5 X  for  x  >  25,  which  corresponds  to  a  distance 
of  about  107  orifice  diameters  in  this  case.  An  acid  solution  of  this  concen¬ 
tration  would  have  a  pH  value  of  about  6,  which  for  most  practical  considera¬ 
tions  is  only  mildly  acidic. 

Class  II  Accident:  Acid  Aerosol  Formation 
The  various  possible  mechanisms  of  acid  aerosol  formation  due  to  an  oleum 
spill  on  water  have  been  reviewed  in  considerable  details  in  the  Interim  report 
(July  1978).  However,  the  actual  scenario  is  too  complicated  to  be  described 
by  a  single  mechanism  or  a  combination  of  several  mechanisms.  The  partition  of 
the  spilled  acid  between  the  bulk  water  and  air  cannot  be  calculated  a  priori. 
Only  some  very  conservative  estimates  could  be  deduced  from  thermodynamic  con¬ 
siderations.  Some  heat  of  formation  data  for  oleum  are  available  in  the  lit¬ 
erature  [23,24].  From  these  data  the  partial  molar  heats  of  mixing  required  in 
the  enthalpy  calculations  may  be  derived.  The  distribution  of  energy  depends 
upon  the  dynamics  of  the  mixing  and  vaporization  processes  and  may  be  expected 
to  vary  considerably  under  different  conditions.  One  can  only  calculate  from 
constant  enthalpy  considerations,  by  assuming  adiabatic  conditions  and  instan¬ 
taneous  and  uniform  mixing,  the  amount  of  SO3  vaporized,  given  a  priori  knowl¬ 
edge  of  temperature  increases  of  the  oleum  which  is  being  diluted  down  to  100% 
H2SO4.  The  results  of  such  calculations  for  20%  oleum  are  given  in  Fig.  13, 
where  the  fraction  of  SO3  vaporized  is  plotted  vs.  the  temperature  Increment. 

It  is  seen  that  as  much  as  12%  of  SO3  could  be  released  into  air  to  form 
sulfuric  acid  aerosols,  if  all  energy  generated  from  dilution  were  channeled 
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into  the  SO3  vaporization  process  only.  Obviously,  this  represents  a  hypothe¬ 
tical  upper  limit;  the  actual  fraction  is  more  likely  to  be  much  less  since 
most  of  the  released  heat  will  be  consumed  heating  up  the  fluid  rather  than 
vaporizing  SO3.  Thus  the  actual  partition  coefficients  have  to  be  measured 
experimentally . 

In  a  typical  oleum  spill  experiment,  the  water  tank  was  initially  filled 
with  tap  water  to  a  predetermined  level  and  the  oleum  release  device  (as  shown 
in  Fig.  1)  was  suspended  from  the  top  center  of  the  tank.  The  tank  was  then 
sealed  against  air  leakage  and  purged  with  nitrogen  gas.  The  release  device  in 
this  case  consisted  of  a  fixed  stainless  steel  spike  and  a  long  spring-loaded 
lucite  plunger  rod.  A  specially  blown  thin-bottom  glass  flask  containing  an 
oleum  sample  of  known  weight  was  capped  and  attached  to  the  lower  end  of  the 
plunger.  When  the  rod  was  released,  the  flask  was  propelled  downward  with  the 
plunger  and  ruptured  against  the  spike  located  a  fixed  distance  from  the  bottom 
of  the  flask.  The  level  of  the  water  in  the  tank  with  respect  to  the  sharp 
point  of  the  spike  determined  the  depth  of  release,  at,  above,  or  below  the 
surface.  Thus,  oleum  was  instantaneously  released  in  a  consistent  manner  at  a 
reproducible  depth. 

Usually,  a  white  cloud  was  observed  to  rise  at  the  instant  of  the 
release.  Four  circulating  fans,  each  located  at  one  corner  of  the  tank  cover 
and  interconnected  through  ducts,  quickly  dispersed  the  cloud  into  a  dense  fog 
which  filled  the  whole  air  space  above  the  water.  Figs.  14a  and  14b  show  the 
fog  formation  in  a  sequence  of  pictures  taken  immediately  before  and  after  an 
oleum  release  experiment.  A  solenoid  valve  controlled  by  a  pressure  switch 
was  then  activated  to  regulate  the  rate  at  which  fresh  water  was  introduced  to 
fill  the  tank.  The  rising  water  level  displaced  a  steady  volume  of  gas  through 
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either  a  membrane  filter  for  acid  aerosol  collection  or  a  diffusion  battery 
(Fig.  15)  for  particle-size  distribution  measurements.  Sulfuric  acid  aerosols 
collected  on  Fluoropore  filters  were  subsequently  analyzed  for  total  sulfate  by 
ion  chromatography  [25].  Fig.  16  shows  the  temperature  variation  with  time  at 
two  local  points  in  water  and  in  air,  indicating  the  extent  of  temperature  rise 
from  oleum-water  reaction. 

Oleum  releases  varying  from  50  to  475  g  in  size  and  containing  either  20% 
or  30%  free  SO3  were  employed.  The  quantities  of  sulfuric  acid  aerosols  pro¬ 
duced  in  these  experiments  were  found  to  depend  on  both  the  sample  size  and 
release  depth  in  water.  In  general,  less  aerosol  was  formed  in  experiments 
where  oleum  was  released  below  the  water  surface  than  at  or  above  the  surface . 
While  the  above-surface  releases  always  generated  a  large  amount  of  both  steam 
and  acid  aerosols,  which  quickly  formed  a  dense  fog  in  the  air  space  of  the 
spill  tank,  the  acid  aerosols  produced  from  the  below-surface  releases  were 
barely  observable  even  with  the  aid  of  a  He-Ne  laser  beam. 

The  parti  "e-size  distributions  were  measured  with  a  "Collimated  Holes 
Structures"  diffusion  battery,  the  operation  details  of  which  have  been  given 
elsewhere  [26,27].  A  nonlinear  iterative  inversion  method  has  been  developed 
[27]  for  the  analysis  of  diffusion  battery  data  (Appendix  5).  Fig.  17  shows  a 
composite  of  two  particle-size  distributions  for  the  acid  aerosols  generated  by 
releasing  50  gm  of  20%  oleum  at  2.5  cm  below  the  surface  (dashed  curve)  and  at 
2.5  cm  above  the  surface  (solid  curve).  Both  sets  of  data,  measured  at  10% 
relative  humidity,  exhibit  a  distinct,  well-defined  log-normal  particle-size 
distribution.  The  particle  sizes  of  the  aerosols  produced  in  underwater 
releases  are  usually  smaller  than  those  produced  in  above-surface  releases. 

A  detailed  study  of  the  effect  of  oleum  release  depth  on  particle  sice  was 
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undertaken  for  20%  oleum  and  the  results  are  shown  In  Fig.  18,  where  the 
droplet  diameter  of  the  acid  aerosol  is  plotted  as  a  function  of  normalized 
sample  immersion  depth.  Here,  an  immersion  depth  is  defined  as  the  distance 
measured  from  the  water  surface  to  the  center  of  gravity  of  the  oleum  sample  at 
the  moment  of  release.  For  the  case  where  the  center  of  gravity  of  the  sample 
is  above  water  surface,  the  Immersion  depth  is  considered  to  be  negative  by 
definition.  The  immersion  depth  is  then  divided  by  the  sample  size  taken  as 
the  radius  of  a  sphere  of  equivalent  volume  to  obtain  the  normalized  sample 
immersion  depth,  Z,  plotted  as  the  abscissa  in  Fig.  18.  A  definite  decrease  in 
particle  size  is  noted  with  increasing  release  depth  until  Z  reaches  about 
2.5.  Below  this  depth,  the  particle  size  seems  to  level  off  at  0.14  i  0.03  urn 
in  diameter.  In  all  cases,  the  observed  aerosol  sizes  (0.1  to  0.6  pm)  are  well 
within  the  respirable  size  range. 

Experiments  were  also  performed  to  determine  the  extent  to  which  sulfuric 
acid  could  become  airborne  in  an  oleum  spill  accident.  Although  such  informa¬ 
tion  is  extremely  useful  in  assessing  the  consequences  of  a  given  accident, 
reproducible  experimental  data  are  very  difficult  to  obtain  since  the  vaporiza¬ 
tion  process  depends  strongly  on  the  manner  in  which  oleum  is  being  released  in 
water.  In  addition,  there  is  no  a  priori  information  to  determine  how  the 
liberated  heat  of  dilution  is  to  be  channeled  into  the  several  processes 
simultaneously  operating  during  the  initial  vaporization  period. 

An  attempt  was  made  to  correlate  the  airborne  acid  fraction  with  release 
depth  and  initial  sample  size.  Fig.  19  shows  the  results  of  two  series  of 
experiments  with  30  Z  oleum  of  two  different  initial  sample  sizes,  namely 
170±10g  and  450±25g,  released  at  the  same  depth,  i.e.,  surface  release.  The 
concentrations  of  the  airborne  sulfate  as  collected  on  filters  are  monitored 
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for  the  duration  of  the  run.  The  concentration  decay  shown  in  Fig.  19  is 
typical  of  a  stirred  aerosol  in  a  confined  space  where  settling  and  plate-out 
are  the  dominant  aerosol  loss  processes.  The  bar  lengths  of  the  data  points 
represent  the  time  duration  that  a  filter  sample  is  taken,  which  ranges  from 
one  to  ten  minutes.  The  percent  sulfate  airborne  for  each  release  was  obtained 
by  extrapolating  the  concentration  to  t“0.  It  is  seen  that  for  the  two  series 
of  experiments,  the  percent  sulfate  airborne  is  almost  Identical  within  the 
precision  of  the  experimental  data,  although  the  initial  sample  sizes  differ  by 
a  factor  of  almost  three. 

Table  2  is  a  summary  of  all  release  data  obtained  for  both  202  and  302 
oleum.  To  normalize  the  release  depth,  the  previously  defined  normalized 
immersion  depth,  Z,  is  again  used  in  Table  2.  Although  data  are  somewhat  scat¬ 
tered,  a  few  conclusions  may  be  drawn  from  the  results.  It  seems  that,  for 
each  type  of  oleum,  there  exists  a  limiting  value  to  the  amount  of  sulfuric 
acid  aerosol  that  may  be  formed  during  a  spill.  The  fractional  vaporization 
does  not  depend  strongly  on  the  initial  sample  size,  as  is  suggested  by  the 
data  shown  in  Fig.  19.  As  a  conservative  estimate,  therefore,  one  can  state 
that  no  more  than  one-tenth  of  a  percent  of  sulfuric  acid  may  become  airborne 
to  form  an  acid  mist  under  most  oleum  spill  conditions.  The  fractional  air¬ 
borne  acid  decreases  very  rapidly  from  surface  spills  (Z  -  0)  to  releases  just 
below  the  water  surface  (Z  -  2). 


Table 2.  Summary  of  airborne  acid  fractions  in  oleum  release  experiments 


Release  Condition 

%  Acid  Airborne 

Release  Position 

Z 

Av.  Sample 
wt.  gm 

20%  Oleum 

30%  Oleum 

S3  cm  above  Surface 

-8 

450 

1.4 -1.6xl0*2 

6.0  -  6.5  x  lO*2 

On  Surface 

D 

175 

450 

BBS 

1.8  -  4.6  x  10*2 
1.1  -  5.1  x  lO*2 

9  cm  below  Surface 

+2 

450 

2.9  -  3.8  x  10*5 

1.2  -  2.5  x  10*4 

17  cm  below  Surface 

+4 

450 

7.1  —  8.4  x  10*® 

1.7  —3.6  x  10*4 

CONCLUSIONS 


Two  aspects  of  sulfuric  acid  spill  accidents  have  been  investigated  in 
laboratory  experiments.  It  is  shown  that  the  convective  mixing  of  concentrated 
sulfuric  acid  with  water  can  adequately  be  described,  for  both  the  instantane¬ 
ous  and  continuous  releases,  by  mathematical  models  taking  into  account  the 
variation  of  the  buoyancy  force  arising  from  changes  in  acid  concentration  and 
liberated  heat  of  dilution.  For  oleum  spills  in  which  acid  aerosol  formation 
is  a  potential  safety  hazard,  a  conservative  estimate  of  less  than  one-tenth  of 
a  percent  is  obtained  to  represent  the  amount  of  airborne  acid  under  most  acci¬ 
dent  conditions.  The  fraction  of  airborne  acid,  however,  decreases  very  rapid¬ 
ly  with  increasing  release  depth.  The  acid  aerosols  exhibit  a  well-defined 
log-normal  particle-size  distribution  with  peak  diameter  varying  from  0.1  to 
0.6  ym  (at  70%  R.H.)  depending  upon  release  depths.  This  is  well  within  the 
respirable  particle  size  range. 
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ACID  SPILL  FACILITY 


Fig.  2  Schematic  diagram  of  BNL  Acid 
Spill  Facility 


Fig.  3  Glass  hemispherical  cups  used  in 
acid  spill  experiments 


Fig.  4  Shadowgraphs  taken  at  various 
stages  of  mixing  for  a  500  g  acid  dump. 
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Fig.  6  Calculated  temperature  and  acid 
concentration  changes  with  time. 


Fig.  7a  Shadowgraphs  of  an  acid  plume 
taken  at  intervals  of  one  frame  per  sec. 


Continued  from  Fig.  7a. 


Fig.  8  Enlargement  of  the 
sec  after  release. 


Fig.  9  Acid  plume  at  21  sec  after  release. 


Fig.  11  Temperature  measurements  at  four 
different  locations  inside  the  acid  plume. 
The  top  insert  was  taken  at  2  inches  below 
the  release  nozzle;  the  bottom  insert  at  9 
inches . 
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Fig.  12  A  comparison  of  the  experimental 
plume  axis  temperatures  (open  circles)  with 
theoretical  prediction  (solid  line).  The 
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Fig.  13  Estimated  fraction  of  SO^ 
vaporized  as  a  function  of  water 
temperature  increase  resulting  from 
heat  of  mixing  (20%  oleum) 


Fig.  14a  Pictures  taken  of  an  oleum  release 
experiment  to  show  the  formation  of  a  dense 


Fig.  15  Picture  of  the  diffusion  battery 
used  In  particle  size  measurements. 
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Fig.  17  Sulfuric  acid  aerosol  size 
distributions  measured  with  the  diffusion 
battery. 
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Fig.  18  Average  droplet  size  of  sulfuric 
acid  aerosols  as  a  function  of  release  depth 
in  water. 


Appendix  1.  Computer  Program  for  Instantaneous  Releases 


The  program  solves  the  dimensionless  governing  equations  for  the  instanta¬ 
neous  release  of  96%  H2SO4  in  stagnant  water.  No  input  data  is  necessary.  The 
output  tabulates  the  following  dimensionless  quantities  from  the  time  of 
release  up  to  a  dimensionless  time  of  600:  time,  volume,  radius,  velocity, 
mass  fraction,  temperature  increase  (°C),  buoyancy,  center  distance,  and 
frontal  distance.  A  description  of  the  FORTRAN  parameters  is  given  as  follows: 
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PROGRAM  01FE 


74/74 _ OPT*  1 


FTN  4,5*414 


PROGRAM  DIFF (OUTPUT) 

DIMENSION  Z{5)  *ZP(5>  .TEMPOS) 

- external- der - 

COMMON  6A,AL1.AL2,AK,S,ETA1,ETA2,ETA3,SHI1,SHI2 
~DATJ  Y07CPT*TB,C70796TI75T63. 091.128, 698*556. 288/ 
AK  ■  0 i 235 


OP  «  1 . 1093A6763-Oi»791  366907*  YO 

ETAl  «  1.0/CP 

ET#2  ■  0.60262078T7CP 

£TA3  ■  1.1093A6763/CP 

- SHI  1  *  6  ,f <5 1 366967/CP - 

SHI2  •  0.1915A9227/CP 


6A  ■  a»yo/CP 
AL1  ■  B*Y0/A 
~ AC2  «X*yo*YO/A 
TIO  «  5.0 

- S  "«1V0  - - - - 

N  «  5 


T  -  0.0 
2(1)  *  1.0 
2(2)  *  0.0 

- ZP)  *  TIO/GA  - -  - 

2(A)  *  1.0 


*  (5 )  *  0.0 

0  *  0.0 

~  PRINT  ”1 00 

100  F0RmaT(1H1,1X,»TIME*,6X.*V0lUME*,9x,*RADIUS*,8X,*vEL0C1TY*,6X, 
♦“  AMASS  FRACT«7AKT*TEMF  INCRrOCT*.AXr*BuOYANCYA,8X,*CENTER*,10X. 
«  «XFRONT«/) 


50  GALL  L'GU(T.H,N,Z,ZP*O.DEK,TEMP) 

52  CONTINUE 

RATT  *  2(l)^*Tl.o/370)  -  "■ 

VEL  o  Z(2)/Z{1) 

- r^^0*Z(4T/2  71) - - - - 

lF(Y.LE.o.tSZ5)GO  TO  20 


TMl  *  EtA3-SHI1*Y 

gopo^o - - - - 

THI  B  ETAl-SHIl«Y 
“GO  TO  AO 
THI  *  ETA2-SHI2*Y 

-continue: - 

T I  ■  GA#z (3) / (ThI*Z ( 1 ) ) 


LUiimmuiM 


PHI  «  1.0* (3.0943-A,083A*ZHI*1.9032*ZHI*ZHI)*ZHI 

-r-e“S“-"070eriT535977*TI - 

AClO  *  PHI *F*Z (4 ) 


XFRONT  *  Z(5)  *  RAD 

- PRIfTTTOl  ,  T  »Z'(1T»RXD7V1 

101  F0RMAT1F6.X.8E15.6.F6.3) 

- 1FTT7GTV6DUT0TGO  TO  51 

DO  10  1*1*5 


7T2T51TXFRONT.ACID 


CONTINUE 
“COP’D  “5  2“ 


I  )  »  ri  9  ™  •  £  f 


lOxunfiJi 


subroutine  oer 


74/74  OPT»l 


FTN  4,5*414 


1 Subroutine  DERmN.z.ZPi 

DIMENSION  2(5) ,ZP(S) 


XOHmON  ~Qr,  AL  f.*L2 ,  AK ,  S ,  ETA  I7ETA27E  TA3 ,  SHI  1 .  SW2 
data  YQ,cp,A,B.c/o,96.i.oi63.o9it 126,698. 558.268/ 


5 

Y  ■  2 /Z ( 1 ) 

VP  ■  V0*V 

IFCYP.LE, 0.6625)80  TO  10 

IF  (YP.LE. 0.8448) GO  TO  20 

lo 

THI  ■  ETA3-SHI 1»YP 

BETA  ■  -SHI1«Y0 

10 

GO^lO— 30 

THI  ■  6TAl-SHIl«YP 

Beta  ■  -shii#yo 

60  TO  30 

15 

20' 

TRI  ■  “ETR2^SHI2«YP 

BETA  ■  -SHI2*Y0 

- 30- 

— CONTINUE 

ZHI  ■  Y/  15 .444 05/ Y 0-4,44* 05* Y) 

20 

PHI  ■  1.0*  <3. 0943-4. 0B34«ANI*1.90'32#2mR2HIT5ZHl 

TEMP  ■  6 A* 2 (3) / (THI*Z ( 1 ) ) 

0  *  OTu  0053597 7*T£MP 

F  ■  S  -  0 

U  «  «  6/F 

P  ■  5.A4405«ZHi*ZHI»t3.0943-B,l66«*ZHI*5.7096*ZHl*ZMI»/<Y0*Y»Y* 

15 

X  PHI ) -  BEtA*a/tHl 

R  ■  -Z(1)*(1.0-AL1*Y*AL2*Y*y)-Z(3)*P 

zp  n )  ■  3 .  o*ak*zT2T7  zn  i  *  *  <r.o  /  3 .0  > 

*1  ■  -H/Z(l» 
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V  i  i 


,,  ill fox' ^’^h’2] 

ax'  ax’  L 

-  f^L  +  y  2r'2  ,  db’l  e"r' 

j_dx'  X2b'  ^  dx'J 


-r'2/X2b'2 


Therefore, 


°o  /  \  /  \  »  r°° 

/2Trr'u,P'(i^i]  (iZ_ )  dr’  -  -irp^u  /  2r'(A-By'+Cy’2)  . 

\3y f  J  f  J 


1  rf  ^ 

dy  ~  2r' 2  db' 1  '(1+T*V*  , 

_L_  +  y  __ — „ e  dr ' 


— —  +  y  -» — T  — -  e 
dx'  X2b' 3  dx' 


In  general, 


fo  fe 


n  r  *  ^ 

1  -(i+“ 

^  ~  2r'2  db’  X7  b77  , 

+  y  __ — e  dr' 

Xzb'  dx'J 


X2b' 2  fdy  ^  2y  1  db’ 
n  dx'  X2+n  b'  dx' 


Xz+n  dx'  Xz+n 


Thus , 


f  2irr'u'p '  (lB)  (®L)  dr'  -  +  s(?)iZ  ££ll 

J  \3y'/  \3x'/  w  L  dx'  b'  dx'J 


where  Q(y)  ”  -4 —  -  —2 —  y  +  — £ —  y2 
X2+l  X2+l  X2+3 


B  ~  C  ~2 

, - -  y  +  — - - -  yz 


(X2+l)2  (X2+2)2  (X^+3)7 
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Appendix  4.  Computer  Program  for  Continuous  Releases 
The  program  solves  the  dimensionless  governing  equations  for  the  continu¬ 
ous  release  of  96%  H2SO4  into  stagnant  water.  No  input  data  is  necessary.  The 
output  tabulates  the  following  dimensionless  quantities  from  the  point  of 
release  up  to  a  dimensionless  distance  of  100:  depth,  flowrate,  radius,  veloc¬ 
ity,  mass  fraction,  temp,  increment.  A  description  of  the  FORTRAN  parameters 
is  given  as  follows: 
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.  #  UAI 


PROGRAM  OIFF 


76/76  OPT=l 


FTN  4.5*414 


1 


5 


15 


25 


25 


3-> 


35 


-G 


45 


5  J 


55 


PROGRAM  DIFF ( OUTPUT) 

DIMENSION  ZUJ,ZPC4>  ,T£MP«4,5» 

PEAL  M3 
EXTERNAL  CfcR 

COMMON  GA,AL1,AL2,  ALPHA,  ETA1,ETA2  ,ETA3, SHU,  SHI  2,  ALSO,  Til 
DATA  G,y37a,  3,C/98i  .u  ,l.96,6  3.C91,128.69a,5'58.2  88/ 

ALPHA  =  .1.11 
AL  =  1.16 
ALSO  =  AL*AL 

CP  =  1. 12  9346 76  3-6,791 3 66  9i!7*T0 

_  ETA1  =  1.4 /CP  _ _ _ _ _ 

"  cYA2  =  i". 60 26267 87/CP" 

ETA3  =  1.1  9346763/CP 
SHI1  =  0.791?669;:7/CP 
SHI2  =  L . 1S1549227/CP 
GA  =  A*V0/CP 

■ _ AL1=_  e*yt/ A _ 

AL2  =  C*Yf*Y./A 
*V  =  20  69.5 
TI.  =  5 .  f 

CX  =  (Mr/G)**  <l.'/3.l'  > 

R  AO  =  i.  ,15/CX 

_ PRINT_  22  (• _ 

ZlC  FORMAT (1H1) 

PRINT  lu?  _  _ 

ICC  FORMATJ//1X,*OEPTH*,5X  ,*FLCW  RATE*,8X,*RADIUS«,7X,*tf£LbCITY*,6X, 
♦  *M  AS  S  FRACT*,4X,*TEMP  INCRC4iC)*1_4X,*BUOYANCY*/» 

N  =  4 

_ h  *  ;_.i _ _ 

zTi>  =  i.t 

_ Z  <2 >  =  l.t _ 

Z  (3 1  =  TIl  /GA 
Z (41  =  l.C 

x  =  1  ,c 

Q  -  * .4 _ 

TI  =  TIC 

_ Y  =  Yf _ _ _ _ _  _ _ __ _ _ 

IFIX.LE.3.l>ALPHA  =  u  *L  £5 
CALL  CEOCX,H,N,Z,ZP,Q,DtR,TEMP) 

101  CONTINUE 

VEL  =  1 ( 2)/ Z(l» _ 

Y  =  YC*Z<4i/Z«i) 

IF<V.  LEj.  4^662  51 GO_TO  cL _ _ 

IFtY.LE.Q.  84481  GO  TO  3<J 
THI  =  ETA3-SHI1*Y 
GO  TO  41 

24  THI  =  £TA1-SHU«Y _ 

GO  TO  40 

31 _ THI  =  _E  TA2  ”  S  H 1 2  *  Y  _ 

4f-  CONTINUE 

TI  =  GA,Z<3I/<THI'*ZC1) > 

50  C  ONTI  NUE 

RAO  =  Z  <i)  /SQRT  <Z  I2H  _ 

p RTNT  1  ,  X, ZTTT ,  R  A  0 , VEL, V, T I 
1  FORMA T(F6.1, 6E15.6) 

I F ( X. GT . 11 i. j) GO  TO  1U2  ,, 

00  13  1=1,5 


program  oiff 


76/76  OPT  = 


FTN  4.54414 


ALPHA  =  t.li  ...  _ 

IFCX.LE.3.0ALPHA  s  2.C65 

_C ALL  D E 0 <  X  , H  ,  N *  Z ,  ZP  ,  Q  *  DE R  * TE MP_I _ 

10  C  ONTI  NUE 

GO  TO  151  - - - - 

1C 2  CALL  EXIT 

_  _ END _ 


“SUBROUTINE  OER  76/76  OPT=l  .  '  FTN  4.54414 


1_  _ SUBROUTINE  OER < X ,N, Z, ZP)  __ 

DIMENSION  Z  (4  )  »  ZP  14) 

_ COMMONJGA, ALL, AL2, ALPHA, ETAl,ETA2,ETA3jSHIl,SHI2, ALSO, TI6 

DATA  G.Yi  , A, B,C/98»,.C  *4.96*63. 091*128. 698, 558. 288/ 

_§ _ T  )/Zfl> _ 

4T  CONTINUE 

_ YP  =  VO  ♦  Y _ 

IF<YP.LE.I .6625IGO  TO  It 

_ IMYP.LE. t . 6448 >GO  TO  2C _ 

10  THI  =  £T A3-SHI1*YP 

_ GO  TO  3i; _ 

__  THj  =  £TA1-SHI1*YP 

_ _ GO  TO  30 _  _  _ _ 

2t'  THI  *  ETA2-SHI2*YP 

15  3£  _  CONTINUE  __  _ 

ZHI  =  Y/<5.444C5/Yl— *.«*44C5*YI 

_ PHI  =  1.,  4<3.;  943-4.  834»ZHI41.9:  32»ZHI*ZHH*ZHI _ 

F  =  l.t  -  t.r,ii,535977*GA*Z(3)/ (THi*7tin 

_ _  5t _ CONTINUE _  _  _ _ _ 

20  U  =  l.a/CALSO  4  l.LI 

v_f.  ali*y/(alsq_4_2.  :» _ _ _  _  _ 

R  =  AL2*Y*Y/<  ALSC  4*3.0 
Q  =  U  -  V  4  R 

S  =  U/ ( ALSO  4  i.C)  -  Y/(ALSQ  4  2.C)  4  R/(ALSQ  4  3.0* 

25  ZP<1)  =  2.'*ALPHA*SQRT(Z<2)> 

ZP 12)  =  2.)  *ALSQ4zm*Za)MPHI4F-i.C  V/2I2>  ~ 

ZPC3)  _*  (  ALSQ+l  •  ZfV*  ((Q-2.{*S)*ZP«l)j4S*Zll)4ZPI2)/ZC2)  ) 

ZP<4)  =  C  .if  *  . * 

RETURN 

30  zNO 

t 

I 

! 


SY  MB  CL IC  REFERENCE  MAP  <R=1> 

POINTS 

OEF 
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Appendix  5.  A  non-linear  Iterative  Algorithm 

for  Inverting  Diffusion  Battery  DATA 
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Abstract 


A  non-linear  iterative  algorithm  la  presented  for  use  in  retrieving 
particle  size  information  from  diffusion  battery  data.  It  is  shown  with 
numerically  generated  test  cases  that  the  inversion  method  has  a  relatively 
fast  rate  of  convergence,  good  size  resolution,  and  is  stable  to  as  much  as 
SZ  random  errors  in  the  data.  Application  of  the  method  to  experimental 
monodi8perse  aerosol  data  gives  satisfactory  agreements  between  retrieved 
particle  sizes  and  corresponding  stated  sizes  up  to  0.8  u®  in  diameter. 
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I.  Introduction 


Diffuslonal  deposition  In  a  flow  stream  through  narrow  channels  has 
long  been  recognized  as  an  effective  means  of  size  classification  of  poly- 
disperse  submlcron  particles.  Since  smaller  particles,  whose  diffusion 
coefficients  are  larger  than  those  of  the  larger  particles,  are  readily 
removed  from  the  stream  by  diffusion  to  the  channel  walls,  diffusion  batter¬ 
ies  consisting  either  of  circular  or  rectangular  channels,  have  been  con¬ 
structed  according  to  this  principle  [1J.  More  compact  and  versatile  ver¬ 
sions  made  of  a  series  of  stainless-steel  screens  are  now  commercially 
available  [2,3]. 

However,  difficulties  in  the  analysis  of  diffusion  battery  data  have, 
until  recently,  prevented  it  from  being  widely  used  for  routine  particle 
sizing.  Several  investigators  [4,5]  have  assumed  a  size  distribution  con¬ 
sisting  of  histograms  over  a  finite  number  of  discrete  size  classes  to  fit 
the  data,  whereas  Raabe  [6,7]  has  adopted  a  lognormal  distribution  to  match 
grouped  data.  Mercer  and  Greene  [8]  have  assumed  instead  that  the  diffusion 
coefficients  of  the  particles  are  lognormally  distributed,  from  which  an 
"equivalent  diffusion"  lognormal  distribution  can  be  deduced  for  the  parti¬ 
cle  diameters  via  an  empirical  correlation  between  diameter  and  diffusi- 
vity.  Recently,  Busigin,  Van  der  Vooren  and  Phillips  [9]  have  assumed  that 
the  actual  size  distribution  can  be  represented  by  a  sum  of  lognormal  dis¬ 
tributions  with  different  parameters  and  adopted  an  iterative  random  search 
technique  to  determine  the  parameters  that  give  a  best  fit  to  the  data. 

In  the  present  paper,  a  new  algorithm,  belonging  to  the  class  of  non¬ 
linear  iterative  inversion  techniques,  is  shown  to  be  well  suited  for  the 
retrieval  of  particle  size  information  from  diffusion  battery  data. 
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AlChough  the  formulation  of  the  algorithm  for  solving  inverse  problems  in 
general  has  been  given  elsewhere  [10],  a  brief  description  of  the  method  as 
applied  to  diffusion  battery  data  is  presented  here  in  section  II.  In  sec¬ 
tion  III,  its  performance  is  checked  with  a  few  test  cases  of  known  size 
distributions.  Finally,  in  section  IV,  it  is  applied  to  actual  experimental 
data  taken  with  latex  spheres  and  sodium  chloride  particles  whose  diameters 
are  known  beforehand. 

II.  Theory 

A  M-stage  diffusion  battery  may  be  depicted  schematically,  as  follows. 


g*(l)  g*(2)  g*(M-l)  g*(M) 


where  f(d)  is  the  unknown  size  distribution  of  an  aerosol  sample,  g*(l), 

1*1 . .  is  the  number  of  particles  which  penetrate  to  the  1th  stage  and 

K(i,d)  is  the  fractional  penetration  of  particles  with  diameter  d  at  the  1th 
stage.  Then,  f  and  g*  are  related  by 

eo 

g*(i)  -  /  K(i,d)f(d)dd  l-l,...,M  (1) 

0 

which  is  the  Fredholm  integral  equation  of  the  first  kind.  To  find  a  solu¬ 
tion  of  Eq.  (1)  is  a  well-known  ill-posed  problem.  In  general,  a  unique 
solution  may  not  even  exist.  There  are  additional  problems  associated  with 
instabilities  in  numerical  computations.  Many  methods  for  obtaining  an 
approximate  solution  have  been  proposed  (see,  for  example,  [11-15]);  each 
has  its  shortcomings  and  advantages  over  others  for  a  specific  class  of  the 
kernel  function,  K(i,d). 
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In  actual  Inversion,  Eq.  (1)  Is  usually  discretized  using  an 
Integration  quadrature, 

N 

g*U)  ’  £ «k*<Mk>*<*k>  1*1,  •  •  •  ,M  (2) 

where  is  the  weight  factor  of  the  integration  quadrature  and  the  diameter 
is  discretized  in  N  points  over  certain  finite  range  [a,b].  The  new 
algorithm  states  that  the  solution  is  obtained  by  iteration  via 

fi(dk)  -  [1+  VlU.dk>!  Vl<dk>  k-l,...,N  (3) 

beginning  with  some  assumed  initial  distribution  f0(dk).  The  corresponding 
penetration  is  then  given  by 

N 

giU)  -  E  WkKU,dk)fi(dk)  1-1 . M  (4) 

k-1 

VlU»dk)  is  chosen  to  have  the  functional  form 

VlU.dk)  -  +  Bi.xCD-Kd.dfc)  +  Ci-xCl)  (5) 

where  A^_i(l),  and  Cj_i(t)  are  coefficients  to  be  determined  by 

some  specified  conditions.  In  order  that  the  Iteration,  Eq.  (3),  always 
produces  positive  values  for  f^(dk),  ¥i~i(t»dk)  is  required  to  satisfy 

|  VlU.dk) |  <  l  (6) 

for  all  1,  t,  and  k. 

It  can  be  readily  proved  from  calculation  that  K(l,d)  is  monotonic 
decreasing  both  in  1  and  d.  Then,  since  f(d)  is  positive  for  all  d,  Eq.  (1) 
implies  that  g*(l)  la  also  nonotonic  decreasing.  Thus,  g*  can  be  so 
labelled  that  g*(l)  and  g*(M)  is  the  largest  and  smallest  values,  respec¬ 
tively.  The  fundamental  feature  of  the  new  algorithm  is  that 
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(7) 


N 

£  wkK(i,dk)fi(dk)  -  gl(i)  -  nlg*o> 

k-1 


N 

£  WkKtM.dk^Uk)  -  gl(M)  -  nMg*(M)  (8) 

k-1 


for  all  1,  where  m  and  hH  are  positive  constants  close  to  1.  The  effect 
of  each  iteration,  under  conditions  (7)  and  (8),  is  to  modify  fj(dk)  such 
that  the  corresponding  gj(l),  1-2, . . . ,M-1,  always  lie  within  the  range 
bounded  by  nlg*(l)  and  qMg*(M).  The  introduction  of  and  nM  is  to 
impart  extra  flexibility  to  the  algorithm  when  dealing  with  experimental 
data  with  random  errors.  Substitution  of  Eq.  (3)  into  Eq.  (7)  and  (8)  gives 


N 

£  wkK(l ,dk) l,dk)f i-l(<lk)  “  0  <9> 

k-1 


N 

£  wkK(M,dk)Tj_j(£,dk)f i-i(dk)  -  0  (10) 

k-1 


for  all  i,  and  a  given  l.  From  Eqs.  (5),  (9),  and  (10),  ’*'i-l( !,dk)  can  be 
written  as 


'ii-l(Mk>  -  Ri-i(i)*#i-i<Mk> 


(ID 


where  $£_i(l,dk) 


Y1-lW  -  jfcfeg]  ♦  «»■*>  - 


Qi-l(t) 


Qi-i<*>  ■  .  .  “*  M 
iC"l  |  ,  N 


Yi-iOO  «  -  - - 

hi-l<M.*>  "  nM8*(.^  hj.iCl.t) 

nig*(i) 

N 

hi-l(j,*)  *  £  WkK(J.dk)K2(i*dk)fi-l(dk) 

k-1 

N 

Pi-i(j,*)  -  £  wkK<:Mk>K(Mk>fi-i(dk) 
k-l 

and  Is  as  yet  undetermined.  Obviously,  the  definition  of 

indicates  that 

|*i_l<A,dk)|  ^  1  (12) 

for  all  i,  i,  k.  Then,  Eq .  (6)  implies  that 

| Ri— 1  ( *•)  |  <  1  (13) 

for  all  i  and  i. 

In  order  that  the  iteration  process  is  convergent,  R^_^(l)  must  be 
chosen  subject  to  condition  (13)  such  that  f^(d)  is  modified  at  each  itera¬ 
tion  in  a  direction  as  to  reduce  the  residual  error  between  g£(i)  and 
g*(l).  However,  because  of  the  ill-posed  nature  of  the  problem,  the  solu¬ 
tion  is  also  required  to  be  as  smooth  as  possible.  Mathematically,  this  is 
expressed  as:  find  R^_^(l)  which  minimizes 

r  -  +  xs(fl)  (14) 

where  6^(1)  Is  the  residual  error  defined  as 
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6,(4)  -  1  - 

g*(t) 


1  N 

-  1  "  pji)  g  •k«l.dk)fi(dk) 

-  «i-i(*>  -  *  Ri-i<l>  (u) 

N 

Zi-lCt)  -Z 

k-1 

S(fj)  is  a  functional  of  f^,  which  la  a  measure  of  the  smoothness  of  f^  and 
A  is  a  Lagrange  multiplier.  This  approach  of  determining  Ri-i(l)  is  analo¬ 
gous  to  that  in  the  method  of  regularization.  Several  measures  of  smooth¬ 
ness  have  been  considered  by  different  investigators,  for  example,  the  vari¬ 
ance,  the  sum  of  the  squares  of  first  differences,  or  the  sum  of  the  squares 
of  second  differences,  etc.  In  the  present  formulation,  the  following  defi¬ 
nition  of  the  measure  of  smoothness  will  be  chosen, 

S(fi)  -/  tf1(d)-fi_1Cd)]2dd 
0 

N 

’  £  «k[*i<*k>  -  fi-lUk>l2 
k-1 

-  D^d)  •  Ri_i(l)  (16) 

N 

where  Di-l(*)  «  £  wkf^i-lCi»dk)^i-lC<*k)]2,  **  should  be  pointed  out 

k“l 

that  this  definition  of  the  measure  of  smoothness  is  only  a  matter  of  con¬ 
venience  to  obtain  an  expression  for  Ri-i(t)»  Thus,  S(fi)  may  be  considered 
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as  proportional  to  the  total  amount  of  adjustment  of  fj_i  at  the 
iteration  and  by  requiring  it  to  stay  within  certain  bounds  has  the  effect 
of  restricting  any  drastic  change  in  f^  during  the  process  of  iteration. 
Substitution  of  Eq.  (15)  and  (16)  into  (14)  and  then  equating  dr/dRi_j(i)  to 
zero  gives 


Ri_l(i) 


Zj-iU) 

g*00 


«i-l<») 


2 

+  XD^U) 


(17) 


It  is  clear  from  Eq.  (17)  that  a  large  value  of  X  reduces  the  amount  of 
change  in  f^-i  at  the  i**1  iteration  and  hence  tends  to  maintain  smoothness 
in  fj,  while  a  small  value  of  X  accelerates  convergence. 

To  initiate  the  iteration,  an  initial  distribution  fQ(dk)  must  be 
chosen  such  that  it  satisfies  Eq.  (7)  and  (8).  Since  the  functional  form  of 
the  sought  distribution  is  in  general  not  known  a  priori,  fQ(dk)  is  first 
assumed  to  be  a  flat  profile,  i.e.. 


^o(^k)  “  C  »  k“l,...,N 


so  that  no  bias  is  aimed  at  certain  sizes.  c  is  evaluated  from  Eq .  (7). 

5  .  ^*(1) 

N 

L  wkK<i.«»k> 

k-l 


Then,  f0(dk)  is  obtained  by  adjusting  f0(dk)  according  to 
fo<0k)  ”  I1  +E’K2(M,dk)  +  F)  fQ(dk) 


(18) 


where  E  and  F  are  determined  by  the  two  conditions  go<l)  -  hig*(l)  and  go(M) 
•  hj4g*(M).  Therefore, 


fo(dk) 


1+ 


ll8*(M)-go(M) 


Lho(M,M)-  h0(l,M), 


K2(M,dk)-J^(I,-*) 
T'lg*(  1) 


C  (19) 


hlg*(l) 


where 


«oW 


N 


»kK(M,dk) 

k-1 


It  should  be  pointed  out  that  the  form  of  f0(dk)  defined  in  Eq.  (18)  is 
deliberately  chosen  to  be  similar  to  'f,i-i(*’,dk) .  If  Eq.  (19)  should  result 
in  negative  values  of  fQ(dk)  for  some  dfc,  | nMg*(M)-g0(M) J  is  successively 
reduced  by  temporarily  adjusting  g*(M)  until  fQ(dk)  >  0  for  all  k  and  then 
Eq.  (19)  is  repeatedly  applied  until  g*(M)  returns  to  its  original  given 
value. 

Twomey  [16]  proposed  a  different  expression  for  'i'i-i( ^,dk) 

Vl<Mk)  *  K<Mk)  \— 1-~  ~  (20) 

|_8i-l(*)  J 

Since  for  any  physical  system  K(l,dk)  can  always  be  normalized  to  be  less 
than  or  equal  to  1,  <1'i-i(i>dk)  in  Eq.  (20)  satisfies  the  condition  (6).  He 
applied  Eq.  (3)  with  (20)  to  determine  the  size  distribution  from  measure¬ 
ments  of  the  number  of  particles  through  nuclepore  filters  and  showed  that 
it  yields  better  results  than  the  constrained  linear  inversion  method,  par¬ 
ticularly  when  the  sought  distribution  is  not  smooth.  Although  the  Twomey* s 
scheme  is  simpler,  the  present  technique,  derived  from  a  general  mathemati¬ 
cal  basis,  is  considered  more  systematic  with  built-in  control  of  speed 
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of  convergence  and  smoothness  of  the  solution.  In  the  next  section,  the  two 
methods  will  be  compared  by  means  of  test  cases  of  known  distributions. 

The  iteration  procedure  of  the  present  algorithm  is  described  as 
follows: 

(1)  compute  fjjCdjj)  from  Eq.  (19) 

(2)  begin  iteration  at  1*2 

(3)  put  X  equal  to  some  chosen  value 

(4)  compute  ^i-i(  *->djt) 

(5)  calculate  Ri-i(l)  from  Eq.  (17) 

(6)  test  if  |Ri_i(i)|  <  1;  if  not,  increase  X  successively  by  one 
order  of  magnitude  until  |k^_^(i)|<  1 

(7)  adjust  fi_i(dit)  by  Eq.  (3) 

(8)  compute  gj(l'),  Jl'*2, . . .  ,M-1 

(9)  increase  i  by  1  and  go  back  to  step  (3)  until  t»M-l 

(10)  calculate  dj^i1),  l'm2 . M-l 

(11)  find  max  1 6 1 ( *. * )  I  ;  if  it  is  less  than  certain  specific 

i'*2, . . . ,M-1 

error  limit  criterion,  stop  iteration,  if  not,  increase  i  by  1  and 
go  back  to  step  (2)  to  repeat  the  iteration  process. 

One  pass  of  the  iteration  procedure  from  step  (2)  to  (10)  will  be  called  one 
cycle  of  Iteration. 

III.  Test  Cases 

Numerical  test  cases  are  generated  by  passing  hypothetical  aerosols 
with  known  distribution  through  the  diffusion  battery.  At  Brookhaven 
National  Laboratory,  experimental  research  is  being  conducted  using  a 
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prototype  diffusion  battery  based  on  the  original  design  of  Sinclair  [1]. 
It  is  made  of  collimated  hole  structures  with  approximately  9750  circular 
channels  in  each  stage  and  has  eleven  stages  and  twelve  sampling  ports 
(H"12).  The  lengths  of  the  various  stages  are  listed  in  Table  1. 

Exact  mathematical  expression  has  been  derived  for  the  fractional 
penetration  K(£,d)  by  Gormley  and  Kennedy  [17]. 


Table  1.  Lengths  of  Diffusion  Battery  Stages 


Stage 

Length  (cm) 

1 

0.3175 

2 

0.3175 

3 

0.6350 

4 

0.9525 

5 

1.2700 

6 

1.9050 

7 

2.5400 

8 

3.8100 

9 

5.0800 

10 

7.6200 

11 

10.1600 

K(i,d) 


K( t-l,d) 


g  9x(l*®n) 


(21) 


where 


9(x,B) 


e  -  ex*/2 


£  Jn(B*2>n 

n-1 


Jn  - 


Jn-1 


Jo  “  1 
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and  f$n  are  roots  of  0(1,0)  -  0,  p  -  is  the  diffusion  parameter,  D  the 
diffusion  coefficient  of  a  particle  of  diameter  d,  L  length  of  the  circular 
channel  between  ports  &  and  £-1  and  q  is  the  volumetric  flow  rate  per 
channel.  However,  the  five  terms  approximation  proposed  by  Thomas  [18]  will 
be  adopted. 


K(£,d)  -  K(i-l,d)  [0.819  e"3,65,J  +  0.097  e~22,3li  + 

0.032  e“57u  +  0.027  e“123M  + 

0.025  e"750*1]  (22) 

£*2,...,M.  This  expression  has  been  verified  experimental  by  Sinclair  et 
al.  [19]  for  particles  smaller  than  0.2  pm  in  diameter.  In  this  paper,  it 
is  assumed  to  be  still  valid  for  particles  up  to  1  pm  in  diameter.  The 
diffusion  coefficient  0  is  calculated  from  Einstein  equation  with  the 
Cunningham's  slip  correction, 


D 


kT 

3v£d 


[’ 


1  +  (1.257  +  0.4 


e~0.55d/oj 


where  k  is  the  Boltzmann  constant,  T  the  absolute  temperature  of  the  carrier 
gas,  C  the  viscosity  of  the  carrier  gas,  o  the  mean  free  path  of  the  gas 
molecules.  It  should  be  noted  that  the  sampling  ports  are  so  labelled  that 
the  first  port  (£"1)  corresponds  to  the  input  of  the  diffusion  battery, 
i .  e . , 


K(1  ,d)  -  1  for  all  d  (23) 

Eq.  (23)  indicates  that  g*(l)  represents  the  total  number  of  particles 
entering  the  diffusion  battery.  Henceforth  in  the  discussion  of  test  cases 
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and  experimental  measurements,  the  diffusion  battery  data  g*(i),  1*1,..., M, 
mill  be  divided  by  g*(l)  so  that  g*(l)  becomes  1;  in  other  words  the  inver¬ 
sion  algorithm  will  yield  a  normalized  size  distribution.  Accordingly, 
can  be  fixed  at  1. 

Since  in  many  actual  data  Inversion,  the  significant  size  range  of  the 
aerosol  sample,  whose  size  distribution  is  sought,  is  often  insufficiently 
known  a  priori,  a  diameter  range  [0.001  urn,  10.0  pm]  over  four  decades  of 
magnitude  will  be  assumed  which  is  considered  to  be  more  than  sufficient  for 
most  applications  of  diffusion  battery.  This  diameter  range  is  then ^ discre¬ 
tized  uniformly  on  a  logarithmic  scale,  i.e., 

t  x  <k-l)/(N~l) 

dk-dj/^n  k«l,...,N  (24) 

where  d^  -  0.001  pm  and  djj  »  10.0  pm.  Simpson's  rule  will  be  adopted  as 
the  integration  quadrature. 

Test  cases  are  generated  by  assuming  a  normalized  size  distribution  and 
calculating  g*(f.),  i»l,...,M  from  Eq.  (2).  In  all  cases,  fifty-one  (N-51) 
diameter  points  and  a  flow  rate  of  2  i/min  are  used  in  computing  K(fc,dk). 
Because  of  the  inherent  non-uniqueness  character  of  the  problem,  certain 
measures  have  to  be  devised  to  terminate  the  iteration  process  at  the  "best’’ 
solution.  However,  the  standard  of  the  "best”  solution  is  somewhat  ambig¬ 
uous,  depending  on  the  quality  of  the  given  data  g*  and  the  desired  accuracy 
of  the  retrieved  f.  In  the  following  discussion  of  test  cases,  the  maximum 

value  of  the  absolute  value  of  6^(t),  t«2 . M-l,  after  each  cycle  of 

iteration,  is  calculated, 
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(25) 


6^  “  max  | 6 i ( A) I . 

1-2,..., M-l 

If  6^  is  leas  than  some  prescribed  limit  or  if  the  number  of  cycles  of 
iteration  is  equal  to  a  certain  prescribed  number,  the  iteration  process  is 
terminated.  Our  experience  with  the  present  algorithm  indicates,  however, 
that  as  the  iteration  process  proceeds,  6^  decreases  and  then  in  some  cases 
it  starts  to  increase  again  at  certain  point  of  the  iteration  process. 

Thus,  two  more  provisions  are  incorporated  in  the  computer  program  to  record 
the  two  "best"  iterated  distributions  in  the  iteration  process  according  to 
these  two  criteria: 

(1)  The  retrieved  distribution  at  the  I^h  cycle  of  iteration  which 
gives  the  smallest  compared  with  that  Riven  by  iterated  distributions  at 
iterations  prior  to  and  after  1$,  that  is 

6  -  min  (26) 

Ifi  i 

(2)  The  retrieved  distribution  at  the  cycle  of  iteration  which 

gives  the  smallest  global  residual  error  in  the  least  square  sense 
compared  with  that  given  by  iterated  distributions  prior  to  and  after  I£, 
that  is 

ei£  -  min  \  (27) 

M-l 

where  -  V  6.2(f)  after  each  cycle  of  iteration. 

1-2 


13 


CTg  -  2.0 

X  -  0.1,  \  -  1.0 

Test  case  3  (Fig.  3)  .  A  lognormal  distribution  with  dg  ■  0.8  pm  and 

og  -  1.5 

A  -  0.1,  ifc  -  1.0 

Test  case  4  (Fig.  4)  A  lognormal  distribution  with  dg  “  0.8  um  and 

og  -  3.0 

A  ■  0.1,  hfi  ■  1.0 
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Teat  case  5  (Fig.  5)  a  bi modal  distribution,  which  is  the  sum  of  two 

lognormal  distributions:  dg«*0.03  pm,  Og  ■  1.2  and 
dg  “  0.15  pm,  Og  “  1.4,  with  the  relative  amount 
of  each  distribution  in  the  ratio  1. 0:5.0, 
respectively. 

A  -  0.1,  Hm  ■  1*0 

Teat  case  6  (Fig.  6)  a  bi modal  distribution,  which  is  a  sum  of  two 

lognormal  distributions:  dg  *  0.09  pm,  Og  >  1.2 
and  dg  -  0.6  pm,  dg  ■  1.4,  with  the  -relative 
amount  of  each  distribution  in  the  ratio  0. 5:5.0, 
respectively. 

A  -  0.1,  \  -  1.0 

Teat  case  7  (Fig.  7)  a  gamma  distribution 

f(dk)  -  up  -  . 

4  dk  / 

where  C  is  the  normalizing  factor 
A  -  0.01,  ^  .  i.o 

Teat  case  8  (Fig.  8)  a  Junge  distribution 

...  .  f  0.0  -  0.001  pm  <  dfc  <  0.05  pm 

"  c*djf  o.05  pm  7  dk  <  10.0  pm 

where  C  is  the  normalizing  factor 

i  0.01  Tj  >  10'3 

1  0.001  Tt  <  10"3 

\  -  1.0 

Teat  caae  9  (Fig.  9)  aame  as  Test  case  2,  except  that  52  random  noise 

was  superimposed  on  g*(t),  t»2,...,M,  by  making 
use  of  a  random  number  generator  in  the  computer 
function  library.  In  this  case,  hg  is  randomly 
fixed  at  within  52  about  1.0.  Since  the  amount  of 
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error  in  g*  is  known  a  priori ,  the  iteration  ia 
terminated  as  soon  as  6^  <  5%. 


X  -  1.0 

Test  case  10  (Fig.  10)  same  as  Test  case  4  except  that  5%  random  noise 

was  superimposed  on  g*(l),  i  -  2 . M 

X  -  1.0 


Table  2.  Comparison  of  performance  between  methods. 


Test  case 

Present 

scheme 

Twomey's 

scheme 

X6 

% 

X« 

1 

10000 

8.3  x  10-5 

10000 

7.6  x  10-3 

2 

694 

1.0  x  10-5 

8.4  x  10_A 

3 

1268 

9.7  x  10-6 

10000 

5.8  x  10“3 

4 

1861 

1.0  x  10-5 

2722 

1.1  x  10-3 

5 

1605 

5.0  x  10-5 

1.7  x  10-3 

6 

8226 

5.0  x  10-5 

10000 

8.0  x  10-* 

7 

33 

4.8  x  10-5 

10000 

9.5  x  10-5 

8 

2431 

5.0  x  10~5 

10000 

9.6  x  10-* 

9 

24 

5.0  x  10~2 

73 

5.0  x  10-2 

10 

23 

5.0  x  IQ-2 

107 

5.0  x  IQ-2 

The  ten  test  cases  Just  presented  were  constructed  from  size  distribu¬ 
tion  models  commonly  used  to  describe  atmospheric  aerosol  particles.  Agree-  t 

ment  between  the  retrieved  distributions  by  the  Inversion  scheme  developed 
in  this  paper  and  the  original  distributions  is  excellent  in  most  cases. 

Particularly  in  the  last  two  test  cases,  the  method  was  found  to  be  stable 

to  5Z  random  error  in  the  given  data.  In  comparison,  the  results  by  ' 

Twomey's  scheme  are  in  general  less  satisfactory  with  respect  to  the  present 

scheme  and  it  has  a  much  slower  rate  of  convergence  to  attain  a  given  degree 

of  accuracy,  although  it  is  simpler  to  implement.  Furthermore,  the  present  | 
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scheme  has  an  added  advantage  in  that  it  has  built-in  control  of  the  rate  of 
convergence  and  the  smoothness  of  the  retrieved  distribution.  If  the  origi¬ 
nal  distribution  is  predominant  in  sizes  smaller  than  about  0.1  pm  in  diam¬ 
eter,  Twomey's  scheme  gives  results  comparable  in  quality  to  that  obtained 
by  the  present  scheme,  though  it  is  still  slower  in  convergence  (see  Test 
cases  2,  5,  and  8).  For  distributions  predominant  in  sizes  greater  about 
0.3  pm  in  diameter,  Twomey's  scheme  is  very  slow  in  convergence  (see  Test 
cases  1,  3,  6,  and  7).  Hitherto,  diffusion  battery  has  been  used  mainly  in 
sampling  particles  smaller  than  0.1  pm.  Now,  with  the  new  inversion  scheme 
developed  in  this  paper,  it  is  possible  to  extend  the  particle  sizing  range 
up  to  0.8  pm  in  diameter,  provided  the  flow  rate  is  not  too  high. 

Although  the  test  cases  above  were  based  on  the  collimated  hole  struc¬ 
tures  diffusion  battery,  it  is  of  interest  to  test  the  performance  of  the 
present  inversion  scheme  for  the  parallel  plate  diffusion  battery  con¬ 
structed  by  Busigin  et  al.  [9]  using  the  hypothetical  gamma  distribution 

.  C  /  4  x  10"4  \ 

T  (■“?“) 

as  the  basis  of  comparison  with  their  technique.  Because  of  the  design  of 
this  diffusion  battery,  the  penetration  at  successive  ports  are  0.965, 

0.861,  0.247,  0.023,  and  0.002  for  a  flow  rate  of  2  i/mln.  In  fact,  the 
penetration  at  the  last  port  is  so  small  that  considerable  computer  time 
was  spent  on  calculating  the  initial  profile  via  Eq.  (19).  Accordingly, 
data  from  the  last  port  was  neglected  and  the  result  of  the  present  inver¬ 
sion  scheme  is  shown  as  the  dashed  curve  in  Fig.  11  with  X  "  0.01  and  a  size 
range  [0.01  pm,  1.0  pm].  The  accuracy  at  6396  cycles  of  iteration  is  6  < 
10“5.  Also  shown  in  Fig.  11  (open  circles)  is  the  result  using  the 


colliaated  hole  structures  diffusion  bsttery  (with  the  first  6  ports  only). 
The  sgreeaent  with  the  original  distribution  is  only  fair  for  the  parallel 
plate  diffusion  battery,  whereas  the  agreement  is  excellent  for  the  colli- 
aated  hole  structures  diffusion  battery.  The  reason  for  the  poorer  perfor- 
aance  of  the  present  inversion  scheme  in  the  case  of  the  parallel-plate  dif¬ 
fusion  battery  may  be  attributed  to  the  fact  that  the  penetration  drops 
sharply  at  the  last  three  ports  and  consequently  there  is  insufficient 
Information  to  yield  adequate  resolution.  While  the  technique  of  Buslgln  et 
al.  gives  slightly  better  agreement,  the  present  method  is  faster  (about  11 
seconds  of  execution  time)  compared  to  1  minute  by  their  method. 

IV.  Interpretation  of  Experimental  Data 

The  success  of  the  inversion  scheme  developed  in  this  paper  is  now  fur¬ 
ther  examined  by  applying  it  to  real  experimental  data.  Monodlsperse  sam¬ 
ples  of  latex  and  sodium  chloride  particles  were  generated  by  an  experimen¬ 
tal  setup  shown  schematically  in  Fig.  (12),  which  has  been  described  else¬ 
where  [20].  Latex  particles  of  sixes  0.312,  0.46,  0.6,  and  0.801  urn  in 
diameter  and  sodium  chloride  particle  of  size  0.26  pm  in  diameter  were 
used.  The  penetrations  of  the  collimated  hole  structures  diffusion  battery 
at  a  flow  rate  of  1.5  i/min  were  measured  by  a  Climet  optical  counter  for 
the  latex  particles  and  by  a  TSI  continuous  flow  CMC  for  the  salt  parti¬ 
cles.  The  normalized  penetrations  (relative  to  port  1,  which  is  the  input 
to  the  diffusion  battery)  are  plotted  in  Fig.  (13).  The  solid  lines  ere  the 
theoretical  penetrations  as  calculated  using  Thomas'  approximate  expression 
(Eq.  [22]).  The  agreement  of  the  experimental  data  with  the  theoretical 
curves  can  be  considered  satisfactory.  Irregularities  in  the  data  may  be 
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attributed  to  several  sources,  the  most  noticeable  of  which  are  fluctuation 
in  the  number  density  of  the  sample  entering  the  diffusion  battery,  unstable 
flow  rate,  and  loss  of  particles  through  sampling  ports  during  measurement. 
Some  justifiable  treatments  were  taken  if  necessary  to  correct  the  data 
before  feeding  them  to  the  inversion  scheme. 

The  retrieved  size  distributions  in  all  cases  are  fairly  broad,  far 
from  being  monodisperse .  This  is  partly  due  to  the  irregularities  in  the 
data  and  partly  due  to  the  non-uniqueness  nature  of  the  data  inversion.  In 
Fig.  14,  the  retrieved  size  distribution  for  0.801  pm  latex  particles  is 
shown  as  the  dashed  curve  over  a  size  range  [0.01  pm,  10.0  pm] .  The  accu¬ 
racy  of  the  inversion,  with  X-0.1,  in  1000  cycles  of  iteration  is  Tj*1.5X. 
Accordingly,  instead  of  taking  the  retrieved  distributions  as  they  are, 
further  useful  information  may  be  extracted  from  them  in  the  following  way. 
Since  the  diffusion  battery  is  an  instrument  to  classify  particles  according 
to  the  magnitude  of  their  diffusion  coefficients,  a  mean  diffusion  coeffi¬ 
cient  calculated  from  the  retrieved  distribution  has  more  physical  signifi¬ 
cance  than  the  size  distribution  Itself,  if  the  particle  sample  is  known  a 
priori  to  be  monodisperse.  From  the  mean  diffusion  coefficient,  a  corres¬ 
ponding  mean  diameter  can  be  deduced.  Since  the  diffusion  coefficient  is 
Inversely  proportional  to  the  diameter,  a  mean  particle  diameter,  d,  may  be 
extracted  from  the  retrieved  distribution  in  the  following  way, 

N  1 

E  wk  —  *  retrieved^k 
k-1  dk 

The  <F  for  the  test  particles  used  in  the  experiment  are  listed  in  Table  3. 
The  agreement  can  be  considered  #air. 
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Table  3.  Comparison 

of  retrieved  particle 

sizes  with  stated  sizes  of  test 

aerosols . 

Stated  Size 

Retrieved 

Test  Aerosol 

(pm) 

Mean  Size  (pm)  Z  Error 

NaCl 

0.26 

0.31 

+19 

Latex 

0.312 

0.33 

46 

Latex 

0.46 

0.44 

-4 

Latex 

0.60 

0.53 

-12 

Latex 

0.801 

0.74 

-8 

V .  Conclusion 

The  new  non-linear  iterative  inversion  scheme  developed  in  the  paper 
has  been  established  to  work  very  well  in  analyzing  data  of  the  collimated 
hole  structures  diffusion  battery.  It  has  a  relatively  fast  rate  of  conver¬ 
gence,  good  resolution,  and  is  stable  to  as  much  as  5Z  random  errors  in  the 
data.  Furthermore,  with  this  inversion  scheme,  the  particle  sizing  capa¬ 
bility  of  the  diffusion  battery  can  be  extended  up  to  0.8  pm  in  diameter  for 
a  not  too  high  flow  rate.  Since  the  commercially  available  screen-type  dif¬ 
fusion  battery  is  usually  calibrated  against  the  original  collimated  hole 
structures  diffusion  battery,  the  new  inversion  scheme  presents  a  convenient 
means  of  analyzing  data  routinely  obtained  from  the  screen-type  diffusion 
battery. 
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Test  case  10. 

Inversion  of  parallel-plant  diffusion  battery  data.  —  ■ ..  assumed 

gamma  distribution;  -  retrieved  distribution  by  present  scheme 

using  parallel-plate  diffusion  battery  data;  °  retrieved 
distribution  by  present  scheme  using  collimated  hole  structures 
diffusion  battery  data. 

Experimental  arrangement  for  generating  samples  and  measuring 
penetration  data  of  Latex  and  sodium  chloride  particles. 

Normalized  penetration  data  for  several  Latex  particles.  - 

theoretical  fractional  penetration. 

Retrieved  distribution  ( - )  by  present  scheme  from  measured 

penetration  data  of  0.801  pm  Latex  particle. 
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